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Abstract 

The performance of muon reconstruction, identification, and triggering in CMS has 
been studied using 40 pb^^ of data collected in pp collisions at a/s = 7 TeV at the LHC 
in 2010. A few benchmark sets of selection criteria covering a wide range of physics 
analysis needs have been examined. For all considered selections, the efficiency to 
reconstruct and identify a muon with a transverse momentum pj larger than a few 
GeV/c is above 95% over the whole region of pseudorapidity covered by the CMS 
muon system, | ?/ 1 < 2.4, while the probability to misidentify a hadron as a muon is 
well below 1%. The efficiency to trigger on single muons with pj above a few GeV/c 
is higher than 90% over the full rj range, and typically substantially better. The overall 
momentum scale is measured to a precision of 0.2% with muons from Z decays. The 
transverse momentum resolution varies from 1% to 6% depending on pseudorapidity 
for muons with pj below 100 GeV/c and, using cosmic rays, it is shown to be better 
than 10% in the central region up to pj = ITeV/c. Observed distributions of all 
quantities are well reproduced by the Monte Carlo simulation. 



Submitted to the Journal of Instrumentation 



© 2013 CERN for the benefit of the CMS Collaboration. CC-BY-3.0 license 



*See Appendix|A]for the list of collaboration members 



1 



1 Introduction 

The primary aim of the Compact Muon Solenoid (CMS) Collaboration is to discover physics 
underlying electro-weak symmetry breaking with the favoured mechanism being the Higgs 
mechanism. Many diverse experimental signatures from other potential new physics should 
also be detectable. In order to cleanly detect these signatures the identification and precise 
energy measurement of muons, electrons, photons and jets over a large energy range and at 
high luminosities is essential. 

In this paper we report on the performance of muon reconstruction, identification, and trig- 
gering evaluated using the data collected by the CMS detector at the Large Hadron Collider 
(LHC) at CERN during 2010. During that period the CMS experiment recorded a sample of 
events produced in proton-proton collisions at a centre-of-mass energy of ^/s = 7 TeV with an 
integrated luminosity of 40 pb^^. Muon reconstruction in CMS has been previously studied 
in great detail using muons from cosmic rays [1, 2J. The first studies using 60 nb^^ of 2010 
proton-proton collision data were reported in Ref. [3|. 
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Figure 1: Longitudinal layout of one quadrant of the CMS detector. The four DT stations in the 
barrel (MB1-MB4, green), the four CSC stations in the endcap (ME1-ME4, blue), and the RFC 
stations (red) are shown. 

A detailed description of the CMS detector can be found in Ref. [4]. A schematic view of the 
detector is shown in Fig.jl] Muon reconstruction is performed using the all-silicon inner tracker 
at the centre of the detector immersed in a 3.8 T solenoidal magnetic field, and with up to four 
stations of gas-ionization muon detectors installed outside the solenoid and sandwiched be- 
tween the layers of the steel return yoke. The inner tracker is composed of a pixel detector and 
a silicon strip tracker, and measures charged-particle trajectories in the pseudorapidity range 
\f]\ < 2.^ The muon system covers the pseudorapidity region \t]\ < 2.4 and performs three 

right-handed coordinate system is used in CMS, with the origin at the nominal collision point, the x axis 
pointing to the centre of the LHC ring, the y axis pointing up (perpendicular to the LHC plane), and the z axis along 
the anticlockwise-beam direction. The pseudorapidity (/ is defined as rj = — lntan(0/2), where cos9 = Pz/p- The 
radius r is the distance from the z axis; the azimuthal angle cp is the angle relative to the positive x axis measured in 
the x-y plane. 
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1 Introduction 



main tasks: triggering on muons, identifying muons, and improving the momentum measure- 
ment and charge determination of high-pj muons. Drift tube (DT) chambers and cathode strip 
chambers (CSC) detect muons in the rj regions of < 1.2 and 0.9 < \rj\ < 2.4, respectively, 
and are complemented by a system of resistive plate chambers (RFC) covering the range of 
< 1.6. The use of these different technologies defines three regions in the detector, referred 
to as barrel (|f/| < 0.9), overlap (0.9 < |f/| < 1.2), and endcap (1.2 < \t]\ < 2.4). Muon energy 
deposits in the electromagnetic calorimeter (ECAL), hadronic calorimeter (HCAL), and outer 
hadronic calorimeter (HO) are also used for muon identification purposes. An event in which 
four muons were reconstructed involving all main CMS subdetectors is shown in Fig.|2] 




(a) (b) 
Figure 2: (a) The longitudinal (r-z) and b) the transverse (r-cp) views of a collision event in which 
four muons were reconstructed. The green (thin) curves in the inner cylinder represent tracks of 
charged particles reconstructed in the inner tracker with transverse momentum pj > 1 GeV/ c; 
those extending to the muon system represent tracks of muons reconstructed using hits in 
both inner tracker and the muon system. Three muons were identified by the DTs and RFCs, 
the fourth one by the CSCs. Short black stubs in the muon system show fitted muon-track 
segments; as the z position is not measured in the outer barrel station, the segments in it are 
drawn at the z centre of the wheel, with their directions perpendicular to the chamber. Short 
red (light) horizontal lines in the r-z view indicate positions of RFC hits; energy depositions in 
the ECAL and HCAL are shown as red (light) and blue (dark) bars, respectively. 

For the measurement of muons the single most important aspect is the choice of magnetic field 
configuration. The requirement of a good momentum resolution, specified to be c(px) / Pt ~1% 
at 100 GeV/c and ~10% at 1 TeV/c, without making stringent demands on spatial resolution and 
the alignment of muon chambers leads to the choice of a high magnetic field. CMS chose a high- 
field solenoid. The magnetic flux generated by the central field is large enough to saturate a 
sufficient amount of steel in the return yoke to allow the installation of four muon stations. This 
provides a good level of redundancy in the measurement of muons. The favourable length to 
radius ratio allows efficient muon measurement up to pseudorapidity of \ri\ < 2.4. The strong 
magnetic field also enables an efficient first-level trigger with an acceptable rate. 

The CMS experiment uses a two-level trigger system. The Level-1 trigger, composed of custom 
hardware processors, selects events of interest using information from the calorimeters and 
muon detectors and reduces the read-out rate from the 20 MHz bunch-crossing frequency to a 
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maximum of 100 kHz O. The high-level trigger (HLT) is software-based and further decreases 
the recorded event rate to around 300 Hz by using the full event information, including that 
from the inner tracker [6J. 

Several types of triggers implemented for the 2010 data taking have been used for the present 
studies. These are discussed in Section |2j together with the resulting data and simulated sam- 
ples. Muon reconstruction and identification algorithms are described in Section |3] The mea- 
sured distributions of various kinematic variables of selected muons are compared with simu- 
lation in Section|4j Section|5]presents muon reconstruction and identification efficiencies for ex- 
clusive samples of prompt muons, kaons, pions, and protons. Section |6] summarizes the results 
on muon momentum scale and resolution for different muon momentum ranges. Backgrounds 
from cosmic rays and beam-halo muons are discussed in Section[7[ Section|8] describes the per- 
formance of different isolation algorithms. Muon trigger performance is discussed in Section|9] 



Section 10 gives a summary of our conclusions. 



2 Data and Monte Carlo Samples 

The data samples used for the muon performance studies reported in this paper were collected 
with the following types of triggers: 

• The zero-bias trigger, defined by the coincidence of signals in two dedicated beam po- 
sition monitors (Beam Position and Timing for LHC experiments, BPTX) in the same 
bunch crossing. Its rate was kept constant at around 20 Hz throughout the year by 
adjusting the prescale factor to compensate for the rising instantaneous luminosity. 
Events collected by this trigger do not suffer from any muon-detection bias at trig- 
ger level and are used to define an inclusive muon sample for the study of muon 
kinematic distributions as discussed in Section H) 

• Single-muon triggers. Muon candidates are reconstructed at the trigger level using 
information from the muon detectors and the inner tracker. Events containing a 
muon candidate with online-reconstructed transverse momentum pj greater than 
a predefined threshold (luminosity-dependent, 15GeV/c or lower in the year 2010) 
are recorded. All muon triggers used in 2010 data taking covered the full muon de- 
tector acceptance corresponding to | f/ 1 < 2.4. These triggers were the main source 
of intermediate-px and high-px muons during 2010 data taking, efficiently selecting, 
e.g., muonic decays of W and Z bosons. Furthermore, to collect cosmic-ray data 
during breaks in LHC operation, a trigger requiring at least two loosely matched 
segments in the bottom half of the barrel muon system or a single segment in the 
endcap muon system was implemented. The triggers from the top half were dis- 
abled to avoid the need for special synchronization. 

• Muon-plus-track triggers. To improve the efficiency of collecting ]/tp events, a special- 
ized high-level trigger was implemented. This trigger selected events in which the 
muon track can be paired with an inner-tracker track of opposite charge yielding an 
invariant mass close to that of the ]/xp. To sample the efficiencies in the whole pj 
region evenly, multiple instances of the trigger were deployed with different thresh- 
olds on the transverse momentum of the inner-tracker track. In addition, another 
set of specific J/i/? triggers was implemented, using only the muon system for the 
reconstruction of one of the two muons. The muon-plus-track triggers were used 
to measure identification and trigger efficiencies for low-px rnuons, as described in 
Sections H] and H] 
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2 Data and Monte Carlo Samples 



• Jet and missing transverse energy (Ej"^^) triggers. Using calorimeter information, jets 
and missing transverse energy are reconstructed online. Triggers with different 
thresholds on jet transverse energy and £^'^* were implemented. These events were 
used to select a sample of muons that was unbiased by the requirements of the muon 
trigger. 

In addition, a loose double-muon trigger requiring two or more muon candidates reconstructed 
online and not applying any additional selection criteria was implemented, taking advantage 
of the relatively low luminosity during 2010 data taking. This trigger selected dimuons in the 
invariant mass region spanning more than three orders of magnitude, from a few hundred 
MeV/c^ to a few hundred GeV/c^, as shown in Fig. [s] The events collected with this trigger 
were used in both the detector commissioning and physics studies. 




Figure 3: Invariant mass spectrum of dimuons in events collected with the loose double-muon 
trigger in 2010. The inset is a zoom of the 8-12GeV/c^ region, showing the three Y(nS) peaks 
clearly resolved owing to a good mass resolution, about lOOMeV/c^ in the entire pseudorapid- 
ity range and 70 MeV/c^ when both muons are within the range \rj\ < 1. 

All collision data samples studied in this paper were filtered by requiring at least one well- 
reconstructed primary vertex to reduce the contamination from non-collision backgrounds. 
Techniques to further suppress the non-collision backgrounds according to the needs of physics 
analysis are discussed in Section [7[ 

To compare the results obtained in data to predictions, a number of simulated samples were 
produced using Monte Carlo (MC) techniques. All MC samples were produced with the 
CTEQ6l [ZI set of parton distribution functions and different event generators were used de- 
pending on the process considered. Samples of tt and QCD multijet events were generated 
using PYTHIA 6 |8| with the Z2 tune [9|, as well as inclusive muon-enriched samples, in which 
only events containing at least one muon with transverse momentum greater than a given 
threshold were selected at generation level. Samples of prompt ]/ip mesons as well as ]/ip 
particles originating from the decays of b hadrons were generated with pythia interfaced to 
EVTGEN [lOJ. Inclusive W and Z samples and non-resonant Drell-Yan events were produced 
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using the POWHEG l|m event generator, interfaced with PYTHIA for the simulation of parton 
showering and hadronization processes. For W+jets and Z+jets samples with a given number 
of jets, the madgraph I^12|I event generator was used, combined with pythia for showering 
and hadronization. 

Generated events were processed through a full GEANT4-based HH [ill detector simulation, 
trigger emulation, and event reconstruction chain. A realistic misalignment scenario based on 
the knowledge of positions of different elements of the inner-tracker and muon systems was 
used to describe the detector geometry. The positions of the tracker modules were evaluated by 
applying the track-based alignment procedure to a sample of 2.2 million cosmic-ray muons and 
3.3 million minimum-bias events collected in 2010. The residual uncertainties in the positions 
of individual tracker modules were measured to be smaller than 6 fim in the pixel detector and 
less than 10 jim in the silicon strip tracker. The procedure used to emulate the remaining mis- 
alignment effects in the simulation closely followed that employed to align the tracker using 
data, and used as input the module displacements determined from cosmic-ray data available 
at the time of the MC sample production. As a result, the simulated geometry of the tracker 
included coherent displacements and rotations of tracker modules that could bias the track re- 
construction without affecting the local alignment precision. These coherent movements were 
estimated not to exceed 200 }im. The alignment precision for muon chambers was estimated by 
comparison with photogrammetry to be about 500 fim for DT chambers and between 300 and 
600 }im (depending on the ring) for CSCs, in the r-(p plane. The misalignment scenario for the 
muon chambers used in the simulation was consistent with this precision. 

Unless stated otherwise, additional proton-proton interactions in the same bunch crossing 
(pile-up) were not simulated because of their relatively small rate in 2010 (2.7 interactions on 
average). Such pile-up is expected to have a negligible effect on the results presented here. 

3 Muon Reconstruction and Identification 

In the standard CMS reconstruction for pp collisions 121 [l5l, tracks are first reconstructed in- 
dependently in the irmer tracker {tracker track) and in the muon system (standalone-muon track). 
Based on these objects, two reconstruction approaches are used: 

• Global Muon reconstruction (outside-in). For each standalone-muon track, a matching 
tracker track is found by comparing parameters of the two tracks propagated onto a 
common surface. A global-muon track is fitted combining hits from the tracker track 
and standalone-muon track, using the Kalman-filter technique ||T6ll . At large trans- 
verse momenta, px ^ 200 GeV/ c, the global-muon fit can improve the momentum 
resolution compared to the tracker-only fit I2ul5j. 

• Tracker Muon reconstruction (inside-out). In this approach, all tracker tracks with pj > 
0.5 GeV/c and total momentum p > 2.5 GeV/c are considered as possible muon can- 
didates and are extrapolated to the muon system taking into accoimt the magnetic 
field, the average expected energy losses, and multiple Coulomb scattering in the 
detector material. If at least one muon segment (i.e., a short track stub made of DT 
or CSC hits) matches the extrapolated track, the corresponding tracker track quali- 
fies as a Tracker Muon. Track-to-segment matching is performed in a local (chamber) 
coordinate system, where local x is the best-measured coordinate (in the r-cp plane) 
and local y is the coordinate orthogonal to it. The extrapolated track and the seg- 
ment are considered to be matched if the distance between them in local x is less 
than 3 cm or if the value of the pull for local x is less than 4, where the pull is defined 



6 



3 Muon Reconstruction and Identification 



as the difference between the position of the matched segment and the position of 
the extrapolated track, divided by their combined uncertainties Q. 

Tracker Muon reconstruction is more efficient than the Global Muon reconstruction at low mo- 
menta, p < 5GeV/c, because it requires only a single muon segment in the muon system, 
whereas Global Muon reconstruction is designed to have high efficiency for muons penetrat- 
ing through more than one muon station and typically requires segments in at least two muon 
stations. 

Owing to the high efficiency of the tracker-track reconstruction |IT7| and the very high efficiency 
of reconstructing segments in the muon system, about 99% of muons produced in pp collisions 
within the geometrical acceptance of the muon system and having sufficiently high momentum 
are reconstructed either as a Global Muon or a Tracker Muon, and very often as both. Candi- 
dates found both by the Global Muon and the Tracker Muon approaches that share the same 
tracker track are merged into a single candidate. Muons reconstructed only as standalone- 
muon tracks have worse momentum resolution and higher admixture of cosmic-ray muons 
than the Global and Tracker Muons and are usually not used in physics analyses. 

The combination of different algorithms provides robust and efficient muon reconstruction. 
Physics analyses can set the desired balance between identification efficiency and purity by 
applying a selection based on various muon identification variables. In this paper we study 
the performance of three basic muon identification algorithms: 

• Soft Muon selection. This selection requires the candidate to be a Tracker Muon, with 
the additional requirement that a muon segment is matched in both x and y coor- 
dinates with the extrapolated tracker track, such that the pull for local x and y is 
less than 3. Segments that form a better match with a different tracker track are not 
considered. These additional requirements are optimized for low pj (< lOGeV/c) 
muons. This selection is used in quarkonia and B-physics analyses in CMS [18J. 

• Tight Muon selection. For this selection, the candidate must be reconstructed outside- 
in as a Global Muon with the x^/d.o.f. of the global-muon track fit less than 10 and 
at least one muon chamber hit included in the global-muon track fit. In addition, 
its corresponding tracker track is required to be matched to muon segments in at 
least two muon stations (this implies that the muon is also reconstructed inside- 
out as a Tracker Muon), use more than 10 inner-tracker hits (including at least one 
pixel hit), and have a transverse impact parameter |dxy| < 2 mm with respect to 
the primary vertex. With this selection, the rate of muons from decays in flight is 
significantly reduced (see Section|4]), at the price of a few percent loss in efficiency for 
prompt muons such as those from W and Z decays (see Section |5]). The Tight Muon 
selection is used in many physics analyses in CMS, in particular in the measurements 
of inclusive W and Z cross sections llT9ll20l . 

• Particle-Flow Muon selection. The CMS particle-flow event reconstruction algorithm ||2T| 
combines information from all CMS subdetectors to identify and reconstruct in- 
dividual particles like electrons, hadrons or muons. For muons, the particle-flow 
approach applies particular selection criteria to the muon candidates reconstructed 
with the Global and Tracker Muon algorithms described above. Depending on the 
environment of the muon (for example, whether it is isolated or not) the selection cri- 
teria are adjusted making use of information from other subdetectors (for example, 
the energy deposition in the calorimeters). In general, the selection is optimized in 
order to identify muons within jets with high efficiency, while maintaining a low rate 
for the misidentification of charged hadrons as muons. The details of the particle- 
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flow muon selection are described in Ref. p2ll . 

The default algorithm for muon momentum assignment in CMS is called the "sigma switch". 
This algorithm chooses from the momentum estimates given by the tracker-only fit and by the 
global fit. The global fit is chosen when both fits yield muon pj above 200 GeV/c and give the 
charge-to-momentum ratios q/p that agree to within Icr^/p of the tracker-only fit; in all other 
cases the tracker-only fit is taken. 

In addition, CMS has developed specialized algorithms for high-px muon reconstruction and 
momentum assignment. As the muon passes through the steel of the magnet return yoke, 
multiple scattering and radiative processes can alter the muon trajectory. While the former is 
not so important for high-momentum muons, the latter can result in large energy losses and can 
also produce electromagnetic showers giving rise to additional hits in the muon chambers. As a 
consequence, the estimate of the muon momentum at the production vertex can be significantly 
different from its true value. Therefore, several different strategies for including information 
from the muon system have been developed and studied using cosmic rays [2J: 

• Tracker-Pliis-First-Muon-Station (TPFMS) fit. This algorithm refits the global-muon 
track ignoring hits in all muon stations except the innermost one containing hits, for 
reduced sensitivity to possible showering deeper in the muon system. 

• The Picky fit. This algorithm again starts with the hit list of the global-muon track, 
but, in chambers appearing to have hits from showers (determined by the hit oc- 
cupancy of the chamber), retains only the hits that, based on a comparison, are 
compatible with the extrapolated trajectory. 

To further improve the resolution at high pj, mainly by reducing the tails of the momentum 
resolution distribution, combinations of the above can be used. In particular, the Tune P al- 
gorithm chooses, on a muon-by-muon basis, between the tracker-only TPFMS, and Picky fits. 
The algorithm starts with the Picky fit, then switches to the tracker-only fit if the goodness of 
fit of the latter is significantly better. Then it compares the goodness of fit of the chosen track 
with that of TPFMS; TPFMS is chosen if it is found to be better. For high-px muons, TPFMS and 
Picky algorithms are selected by Tune P in most of the cases, in approximately equal amounts, 
while the tracker-only fit is selected only in a few percent of events. For most analyses of the 
2010 LHC data involving high-px muons. Tune P was used for the determination of the muon 
momentum. 

4 General Comparisons between Data and Simulation 

In this section we present data-to-simulation comparisons for two samples of muons: 1) a fully 
inclusive sample of low-px muons collected with the zero-bias trigger, and 2) an inclusive sam- 
ple of intermediate- and high-px muons collected with the single-muon trigger requiring a 
minimum transverse momentum of 15GeV/c. 

Events collected with the zero-bias trigger were required to contain at least one reconstructed 
primary vertex within 24 cm of the geometric centre of the detector along the beamline and 
within a transverse distance from the beam axis of less than 2 cm. The efficiency of this re- 
quirement for simulated pp collisions having at least one reconstructed muon was found to 
be 99%. About 14 million minimum-bias events were thus selected from a total data sample of 
events corresponding to 0.47 nb^^ of integrated luminosity; the contamination from cosmic-ray 
muons in the sample was estimated to be negligible. The corresponding MC sample consists 
of about 36 million minimum-bias events generated using PYTHIA. 
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The sample of events collected with the single-muon trigger with the px threshold of 15 GeV/c 
consists of about 20 million events corresponding to an integrated luminosity of 31 pb^^. Monte 
Carlo samples used for the comparison correspond to about 10 times larger integrated luminos- 
ity and include the simulation of QCD processes, quarkonia production, electroweak processes 
such as W and Z boson production, non-resonant Drell-Yan processes, and top-pair produc- 
tion. The total cross sections for W and Z production were rescaled to match the next-to-next- 
to-leading-order (NNLO) calculations; the cross sections for b-hadron and tt production were 
rescaled to the next-to-leading-order (NLO) calculations (see Section |42] |. 

In both cases the simulation was normalized according to the integrated luminosity of the data 
sample. The uncertainty in the absolute value of luminosity was estimated to be 4% p3|. 

4.1 Classification of muon sources in simulation 

In the range of pj < 30GeV/c, the most abundant source of muons is semileptonic decays 
of heavy-flavour hadrons. This contribution is accompanied by a high rate of muon candi- 
dates arising from light-flavour hadron decays and hadron showers not fully contained in the 
calorimeters. The relative weights of these background contributions are quite sensitive to the 
details of the muon selection. Muons from decays of W and Z bosons dominate the pj spectrum 
in the region pj > 30GeV/c. 

In the simulation, for each reconstructed muon, the hits in the muon system can be associated 
unambiguously with the simulated particle that produced them. This allows the classification 
of reconstructed muons into the following categories: 

• Prompt muons. Here the majority of muon chamber hits associated with the recon- 
structed muon candidate were produced by a muon, arising either from decays of 
W, Z, and promptly produced quarkonia states, or other sources such as Drell-Yan 
processes or top quark production. These individual sources are shown separately 
where appropriate. 

• Muons from heavy flavour. Here the majority of muon chamber hits of the muon can- 
didate were again produced by a muon, but the muon's parent particle was a beauty 
or charmed hadron, or a T lepton. This class of events has been split according to the 
heaviest flavour generated in the event. Hence, beauty includes muons from direct 
b-hadron decays, from cascade b — > c hadron decays, as well as cascade decays of t 
leptons from b hadrons. 

• Muons from light flavour. In this category, the majority of muon chamber hits of the 
muon candidate were produced by a muon arising from a decay in flight of light 
hadrons (tt and K) or, less frequently, from the decay of particles produced in nuclear 
interactions in the detector material. This category includes hadrons whose tracks 
reconstructed in the tracker were mistakenly matched to the muon chamber hits. 

• Hadron punch-through. Here the majority of muon chamber hits of the misidentified 
muon candidate were produced by a particle that was not a muon. "Punch-through" 
(i.e., hadron shower remnants penetrating through the calorimeters and reaching 
the muon system) is the most common source of these candidates, although "sail- 
through" (i.e., particles not undergoing nuclear interactions upstream of the muon 
system) is present as well. 

• Duplicate. If one simulated particle gives rise to more than one reconstructed muon 
candidate, that with the largest number of matched hits is assigned to one of the 
above categories, and any others are labeled as "duplicate". Duplicate candidates 



4.2 Kinematic distributions of muons 
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can arise either from failures of the pattern recognition of the reconstruction soft- 
ware, or from patterns that mimic multiple candidates. 

4.2 Kinematic distributions of muons 

From the 2010 data sample of zero-bias events with a well-established primary vertex, we ob- 
tain 318 713 muon candidates passing the Soft Muon selection and 24334 passing the Tight 
Muon selection. The overall ratio of the number of muon candidates in data to the predic- 
tion of the PYTHIA MC generator normalized to the same integrated luminosity is 1.05 for Soft 
Muons and 1.01 for Tight Muons. 

The distributions of the muon transverse momentum pj multiplied by its charge q, pseudo- 
rapidity ri, and azimuthal angle (p for Soft and Tight Muons in zero-bias events are shown 
in Fig. |4] The pseudorapidity distribution is peaked in the forward region because there the 
minimum pj required to reach the muon stations is lower than in the barrel: in the endcaps 
the threshold in pj is about 0.5 GeV/c, while in the barrel it is about 3-4GeV/c. Overall, there is 
good agreement between data and simulation both in the number of events and in the shapes of 
the distributions. Some discrepancies result from imperfect simulation of local detector condi- 
tions, affecting for example the muon identification efficiency at low pj, as shown in Section |5] 
Furthermore, the leading-order QCD predictions by PYTHIA have large uncertainties. 

Table [ijlists the sources of muons according to simulation. The majority of reconstructed muon 
candidates originate from decays in flight of pions and kaons ("light flavour"). This is partic- 
ularly evident for Soft Muons, while Tight Muons have larger heavy-flavour components. For 
both selections the contribution of muons from heavy-flavour decays increases with pj. The 
Tight Muon selection reduces the hadron punch-through contribution to 0.2% while it is about 
5% in Soft Muons. The measurements of muon misidentification probabilities presented in Sec- 
tion 5.3 confirm that the simulation correctly estimates the probability for light hadrons to be 



misidentified as muons. 

Table 1: Composition by source of the low-px muon candidates reconstructed in zero-bias 
events, according to simulation for the Soft and Tight Muon selections. 



Muon source 


Soft Muons [%] 


Tight Muons [%] 


beauty 


4.4 


22.2 


charm 


8.3 


21.9 


light flavour 


79.0 


55.7 


hadron punch-through 


5.4 


0.2 


duplicate 


2.9 


<0.01 


prompt 


<o.i 


<o.i 



Among all single-muon triggers used in 2010, the trigger with a pj threshold of 15 GeV/c was 
the lowest-threshold unprescaled trigger during a period when most of the data, corresponding 
to an integrated luminosity of about 31 pb^^, were collected. The kinematic distributions of 
muons collected with this trigger have been compared to the Monte Carlo expectations after 
applying a selection on the reconstructed muon pj of 20 GeV/ c, for which the trigger efficiency 
has reached the plateau. The Tight Muon selection applied in this kinematic range has a high 
efficiency for prompt muons, removing most of the background from light-hadron decays and 
hadron punch-through. After the Tight Muon selection, 824 007 muon candidates remain. 

For comparison with these events, the beauty production cross section given by PYTHIA has 
been rescaled to the NLO QCD predictions 124. .25,1 , which were shown to describe recent CMS 
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4 General Comparisons between Data and Simulation 




(j) [rad] (|) [rad] 



Figure 4: Distributions of kinematic variables for a sample of muons selected by the zero-bias 
trigger, for data (points) and for simulation subdivided into the different categories of con- 
tributing muons (histograms). The kinematic variables are the muon transverse momentum 
multiplied by the charge (top), pseudorapidity (middle), and azimuthal angle (bottom). For 
each variable, the left plot shows the distribution for Soft Muons, and the right plot that for 
Tight Muons. The first (last) bin in the cj x pj distributions includes the underflow (the over- 
flow). The error bars indicate the statistical uncertainty, for both data and MC samples. 
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measurements well ||26l |27|. Without this rescaling, an excess in the predicted beauty com- 
ponent is observed in the inclusive px distribution for muon pj lower than about 40 GeV/ c, 
and also in the muon impact-parameter distribution for transverse distances characteristic of 
b-hadron decays. No corrections accounting for differences between the measured and ex- 
pected muon trigger and identification efficiencies have been applied. Such corrections could 
lead to effects of up to 5%, dependent on pseudorapidity, as discussed in Sections |5] and |9j The 
overall ratio of the muon yield in data to the Monte Carlo predictions normalized to the same 
integrated luminosity is found to be 1.02. 
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Figure 5: Distributions of transverse momentum (top left) and pseudorapidity (top right) for 
Tight Muons with pj > 20 GeV/c, comparing data (points with error bars) to Monte Carlo sim- 
ulation broken down into its different components. The last bin in the pj distribution includes 
the overflow. Dips in the f] distribution are due to inefficiencies related to the muon detector 
geometry. The corresponding ratios of data and MC distributions are shown in the bottom row. 
The error bars include statistical uncertainties only. 

The muon transverse momentum and pseudorapidity distributions for pj > 20 GeV/c are com- 
pared to the expectations from the Monte Carlo simulation in Fig. |5] The estimated composi- 
tion of the sample obtained from the simulation-based studies is shown in Table |2] Muons 
from light-hadron decays are predicted to contribute less than 10%, while the hadron punch- 
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through is suppressed to about 1%. The beauty contribution dominates up to muon transverse 
momentum of about 30GeV/c, where the W contribution starts to prevail, leading to a shoul- 
der in the falling pj spectrum. The inclusive muon yield agrees with the expectations within 
a few percent up to a transverse momentum of 50GeV/c. At higher momenta the leading pro- 
cesses are W and Z production, occasionally associated with hard jets. In this pj region, the 
data agree with the predictions within 10%. This has been verified to be fully consistent with 
theoretical uncertainties related to missing higher-order QCD contributions, by comparing the 
MADGRAPH generator used to simulate W and Z with other Monte Carlo programs for the 
W(Z)+jets processes. In conclusion, given the known experimental and theoretical uncertain- 
ties, the agreement between the data and simulation is satisfactory over the entire momentum 
range of pj < 200 GeV/c. 

Table 2: Composition by source of Tight Muons with pj > 20 GeV/ c according to simulation. 



Muon source 


Tight Muons with pj > 20GeV/c [%] 


W (+ jets) 


20.8 


Z/Drell-Yan (+ jets) 


4.7 


top 


0.1 


quarkonia 


0.7 


beauty 


47.6 


charm 


17.4 


light flavour 


7.8 


hadron punch-through 


0.9 


duplicate 


<0.01 



4.3 Muon identification variables 

The basic selections discussed in Section |3] can be further refined for specific purposes using 
additional information available for each reconstructed muon. 



X10 CMS vs ^ 7 TbV, L ^ 0.47 nb " ' ^ CMS ys = 71eV L = 47 nb -' 




# matched stations |d^yl w.r.t PV [cm] 

(a) (b) 

Figure 6: Comparison of data and simulation for distributions of Soft Muons in zero-bias 
events: (a) number of muon stations with matched segments; (b) transverse impact parameter 
dxy of the muon with respect to the primary vertex (PV). The MC distributions are normalized 
to the integrated luminosity of the data sample. The error bars indicate the statistical uncer- 
tainty. 
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The Tracker Muon segment matching is a powerful tool to reject hadron punch-through. The 
number of muon stations with matched segments (see Section |3| is shown in Fig. |6]^a) for Soft 
Muons in zero-bias events. The contamination from hadron punch-through is evident in the 
expectation for one matched station. The probability for a punch-through to be identified as a 
Soft Muon is drastically reduced by requiring matched segments in at least two stations. The 
contribution from low-momentum muons from light-quark decays is also suppressed by this 
requirement. Figure |6]^b) shows the distribution of the transverse impact parameter, where the 
long tail is dominated by pion and kaon decays in flight. Requirements on both variables are 
used in the Tight Muon selection designed to select prompt muons such as those from W and 
Z decays (see Section |3|. 

A few other examples of identification variables are shown in Fig. [7| for Tight Muons with 
Pt > 20GeV/c collected with the single-muon trigger. The transverse impact parameter and 
its significance in Figs. [7|ja) and (b) are useful to select either prompt muons or, by inverting 
the requirement, muons from heavy-flavour decays. The x^/d.o.f. of the tracker-track fit is 
also a good discriminant to suppress muons from decays in flight, as can be seen from the 
composition of the tail of the distribution in Fig. |7|c). The muon isolation is a simple quantity 
to select prompt muons with high purity. The scalar sum of the transverse momenta of tracks in 
the inner tracker and the transverse energies in calorimeter cells (both in the ECAL and HCAL) 
within a cone of radius AR = ^JJK(p)^^^\^^^KrJY = 0.3 centred on the direction vector of the 
muon candidate is calculated, excluding the contribution from the candidate itself. The relative 
combined isolation JJoJ^^, (further discussed in Section jsj is defined as the ratio of this scalar 
sum to the transverse momentum of the muon candidate. Figure |7]^d) shows the distribution 
of the J^o^jnb variable in events with a single reconstructed primary vertex, compared to the 
simulation with no pile-up effects. Isolated muons promptly produced in decays of W and 
Z bosons dominate the region 1^^^^:^ < 0.1. Overall, the agreement between data and Monte 
Carlo predictions for the muon identification variables is good both for zero-bias events and 
for events recorded with the single-muon trigger. 

The accuracy of the propagation of the tracker tracks to the muon system and the perfor- 
mance of the track-to-segment match have been further studied using Tracker Muons with 
Pt > 20GeV/c, with tight selection requirements on the tracker variables only, to avoid possi- 
ble biases. To further purify the muon sample, an isolation requirement I^gj^b < 0.1 has been 
applied. 

The distribution of distance in local x between the position of the extrapolated tracker track 
and the position of the muon segment has been compared between data and simulation, for 
successful track-to-segment matches (distance less than 3 cm or pull less than 4, see Section |3]|. 
The RMS width of residuals is shown in Fig. |8]as a function of the muon-station number, for the 
DT and CSC systems. As expected, the width of the distributions increases with the amount 
of material upstream of the muon station and with the distance over which the track is extrap- 
olated, from the innermost to the outermost muon stations (from MBl to MB4 and from MEl 
to ME4 in the DT and the CSC systems, respectively). The general trend is well reproduced by 
the Monte Carlo simulation, although the increase of the width is a bit larger in the simulation. 
As the outer ring of ME4 is only partially instrumented with chambers (see Fig. [T]), the residual 
for ME4 is not directly comparable with the measurements in the other stations because muons 
traversing the installed ME4 chambers have a higher average momentum. 

The residuals for track-to-segment matches have also been studied as a function of the muon 
momentum and pseudorapidity. The first muon station is the most important in the global 
track reconstruction: it is where the track's sagitta, determined by the magnetic field inside 
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Figure 7: Identification variables for Tight Muons with pj > 20 GeV/ c: (a) transverse impact 
parameter d^y with respect to the primary vertex (PV); (b) significance of the transverse impact 
parameter; (c) x^/ d.o.f . of the fit of the track in the inner tracker; (d) relative combined isolation 
(tracker+calorimeters), with a cone size AR = 0.3, for events with a single reconstructed PV. 
The MC distributions are normalized to the integrated luminosity of the data sample in (a), (b), 
and (c), and to the number of events in the data sample in (d). The last bin in (b) and in (d) 
includes the overflow. Error bars indicate statistical uncertainties. 



the solenoid, is largest, and its measurements are the least affected by multiple-scattering and 
showering effects because the material upstream of the first station corresponds only to the 
inner detectors, the calorimeters, and the magnet cryostat, whereas the stations downstream 
are also preceded by the steel sections of the magnet return yoke. The average material thick- 
ness traversed by a muon reaching the first station is quite different depending on the angular 
region IITSl . In addition, the width of the position residual increases linearly with the propa- 
gation distance of the muon trajectory from the interaction region to the first muon chamber. 
Figure |9]^a) shows the RMS width of the residual of the local x position as a function of muon 
pseudorapidity. A selection on the minimum momentum p > 90GeV/c has been applied to 
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Figure 8: RMS width of residuals of the local x position, given by the position of the muon 
segment with respect to the extrapolated tracker track, as a function of the muon station, for DT 
chambers in the barrel region (left) and CSCs in the endcap regions (right). Data are compared 
with MC expectations. 

remove the bias induced by the trigger threshold (px > 15GeV/c) in the endcap regions. As 
expected from the rj dependence of the distance between the inner tracker and the first muon 
station and of the material thickness, the residual width reaches the maximum values in the 
overlap region, 0.9 < \f]\ < 1.2. Figures |9]^b),|9]^c), and |9]ld) show the momentum dependence 
of the RMS width, separately for the barrel, overlap, and endcap regions. The width decreases 
with increasing momentum because of smaller multiple-scattering effects. The shapes of the 
distributions are well reproduced by the Monte Carlo simulation, although the simulation pre- 
dicts a somewhat larger width. The only exception is the highest momentum bin in the endcap 
region (Fig.|9]^d)), where the distribution of position resolution in data has larger non-Gaussian 
tails than predicted by simulation. This also leads to the discrepancy observed at the extreme 
t] bins in Fig.|9ja). 

We have also examined the distributions of pulls of the local positions and directions in both 
DT and CSC systems. The widths of the pulls were found to be close to unity and no large 
biases were observed, thus demonstrating that the propagation works as expected and that 
the uncertainties are well estimated. The widths of the pull distributions in the simulation are 
about 10% larger than in data. As demonstrated in the next section, such agreement between 
the expected and the measured residuals and pulls is sufficient to obtain a good description of 
the muon reconstruction and identification efficiencies by the simulation. 

High-momentum muons can give rise to electromagnetic showers in the muon system. These 
may produce extended clusters of hits in the muon chambers, which can degrade the quality 
of muon track reconstruction. Hence an accurate simulation and reliable identification of such 
showers are needed. 

We have considered several variables that might indicate the presence of electromagnetic show- 
ers in the muon system, and examined how well they are reproduced by the simulation. Fig- 
ure [To] shows the distributions of two of them: the number of hits reconstructed in the DT 
chamber crossed by a track but not used in the track fit, and the transverse size of the cluster of 
hits around a track. The cluster of hits is defined iteratively, starting from the impact point of 
the extrapolated muon track and successively adding any hit if it lies within < 0.05 rad of 
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Figure 9: RMS width of residuals of the local x position for the track-to-segment match in the 
first muon station (a) as a function of the muon pseudorapidity, with a requirement on mo- 
mentum p > 90 GeV/ c, and (b)-(d) as a function of the muon momentum in different angular 
regions: (b) \r]\ < 0.9; (c) 0.9 < |//| < 1.2; (d) 1.2 < |//| < 2.4. Data are compared with MC 
expectations. 

the hit in the existing cluster with largest radial distance from the impact point. The transverse 
size of a cluster is defined as the maximum distance in the local x-y plane between the impact 
point of the track and any hit in the cluster. The plots are made for two samples of high-energy 
muons (with p reconstructed by the tracker-only fit above 150 GeV/c): 1) collision muons, and 2) 
cosmic-ray muons selected from collision data samples, with a topology similar to that of colli- 
sion muons. Collision muons were selected by requiring that an event has at least one primary 
vertex reconstructed close to the nominal beam-spot position, and that the muon passes the 
Tight Muon selection with additional isolation and track-quality selections to reduce contami- 
nation from hadron punch-through. Cosmic-ray muons were selected by requiring events with 
at most two tracks reconstructed in the inner tracker, one of which was also reconstructed as a 
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Figure 10: Comparison of data and simulation for variables characterizing electromagnetic 
showers, for muons with p > 150 GeV/ c in the barrel region: the number of reconstructed muon 
hits not used in a track fit (left); the transverse size of the cluster of hits around a track (right). 
All distributions are normalized to the number of muons per muon station in the collision data 
sample. 



Global Muon. The momentum spectra of selected coUision-muon and cosmic-muon data sam- 
ples are quite similar, but the cosmic-muon sample provides a larger number of high-energy 
muons. Monte Carlo samples of Drell-Yan dimuon events passing event selection criteria iden- 
tical to those applied to the collision muons were used for comparisons with the data. The vast 
majority of selected cosmic-ray muons are contained in the barrel region, so Fig. 10 shows dis- 
tributions for the DT chambers only. We observe good general agreement between the data 
and MC simulation. 



5 Muon Reconstruction and Identification Efficiency 

The previous section focused on the comparison between data and simulation for various dis- 
tributions of inclusive samples of reconstructed muon candidates. In this section we study 
exclusive samples of prompt muons, pions, kaons, and protons in data to determine the prob- 
ability that such a particle is reconstructed and identified as a muon. 

Throughout this paper, efficiencies are defined in a relative manner, such that the total effi- 
ciency for the entire muon triggering, reconstruction, and identification chain can be calculated 
as the product of the following individual factors: 



— ^track ■ £rec+id " ^iso ' ^trig • (1) 

The efficiency to reconstruct a muon in the inner tracker etrack was measured separately IIT71IT9) 
and found to be 99% or higher within the whole tracker acceptance, in good agreement with 
the expectation from simulations. Given the existence of a tracker track, the combined muon 
reconstruction and identification efficiencies of the different selection algorithms Crec+id can 
be measured using "tag-and-probe" techniques. In this section the tag-and-probe method is 
described, and measurements of the ^rec+id efficiencies are presented. In Section|8| isolation ef- 
ficiencies eiso are calculated from a sample of identified muons. In Section|9| trigger efficiencies 
etrig are defined relative to muons identified offline and, unless otherwise mentioned, passing 
isolation criteria. 
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5.1 Muon efficiency using the tag-and-probe method on dimuon resonances 
5.1.1 Method 

We evaluate the efficiencies for prompt muons by applying a tag-and-probe technique to mu- 
ons from ]/rp and Z decays. Using this technique it is possible to obtain almost unbiased es- 
timates of the efficiencies of the different stages of muon trigger and offline reconstruction. 
Events are selected with strict selection requirements on one muon (the "tag" muon) and with 
a more relaxed selection on the other muon (the "probe" muon), such that the selection applied 
to the probe muon does not bias the efficiency that one wants to measure. The fraction of probe 
muons that passes the selection under study gives an estimate of its efficiency. 

In this section, muon efficiencies erec+id are measured with this technique. The probes are tracks 
reconstructed using only the inner tracker, so there is no bias from the muon subdetectors. 

In the case of the J/t/^ events, combinatorial backgrounds from other tracks in the event are 
generally high, particularly at low pj. An effective way to suppress this background is to 
require that the candidate probe muon has the signature of a minimum-ionizing particle (MIP) 
in the calorimeters. In this way the background can be reduced by about a factor of three 
without using any information from the muon system. The residual background in both ]/xp 
and Z events is subtracted by performing a simultaneous fit to the invariant mass spectra for 
passing and failing probes with identical signal shape and appropriate background shapes; the 
efficiency is then computed from the normalizations of the signal shapes in the two spectra. 

The uncertainty on the fitted efficiency is determined from the likelihood function. As normal- 
izations of signal and background, efficiency of the background, and parameters controlling the 
shapes of the signal and background are all parameters of the fit, the uncertainty includes the 
contributions from the background subtraction procedure. When the background is not negli- 
gible, as in the case of]/xp, the uncertainty on the efficiency obtained by the fit is dominated by 
these contributions. 

For the Z resonance, an unbiased sample of dimuon pairs can be collected efficiently using 
high-px single-muon triggers. For J/i/?, specialized high-level triggers were implemented, as 
described in Section |2] The muon-plus-track trigger used for the ]/ip case does not bias the 
efficiencies related to the muon system, but introduces a small positive bias in the efficiency 
for the Tight Muon selection, which includes quality requirements on the muon tracker track. 
To measure this bias, the efficiencies for these quality requirements alone are extracted using 
another special dimuon trigger that uses only the muon system for the reconstruction of one of 
the two muons. The bias in efficiency, measured to be (0.7 ± 0.1)% in the barrel and (0.3 ± 0.2)% 
in the endcaps, is well reproduced by the simulation and cancels out in the ratio of efficiencies 
from data and from simulation. 

Under certain kinematic configurations muons from J/t/j decays can be close to each other in the 
muon system. This can result in inefficiencies for some muon identification algorithms. To ob- 
tain an unbiased measurement of single-muon efficiencies, a separation criterion is applied to 
the tag-probe pairs: the extrapolated impact points of the two muon tracks on the surface of the 
first muon station must have an angular separation AR = ^ (A//)^ + (A^)^ > 0.5. The impact 
of this requirement on Z — > events is small: as the large opening angle at the produc- 

tion is preserved by the smaller bending in the magnetic field (~0.1 rad for a px of 25GeV/c), 
only 0.2% of the Z — > events fail the separation criterion above. A dedicated measure- 

ment of the dimuon efficiencies as a function of the separation between muons is described in 
Section l52l 



5.1 Muon efficiency using tine tag-and-probe metlnod on dimuon resonances 
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5.1 .2 Results 

Figure ITT] shows the muon efficiency ej-ec+id given that a tracker track exists, measured using 
]/ip -^jrji^ and Z — > ji^ji^ events. The results obtained from the data collected in the 2010 
LHC data-taking period are compared with those from simulated events. 

For comparisons with Z — )■ events, an unweighted sample of simulated events corre- 

sponding to an integrated luminosity of ~330pb^^ is used: the simulated samples are Z — > 
W+jets, and muon-enriched QCD (see Section|2|. For studies at the ]/ip peak, separate 
samples of prompt ]/xp —> and B —> ]/ip + X —> ji^ji^ + X are used, simulated as de- 

scribed in Section |2] All MC samples used for the results in this section included simulation of 
pile-up. Simulation of the background processes is not included for the ]/tp case, as it would be 
impractical to simulate a sufficient number of inclusive muon-plus-track events. For studies of 
systematic uncertainties described below, samples of background events have been generated 
according to the background invariant mass spectra determined from fits to the ]/xp — > ji^ji^ 
events in the data, and added to the simulated signal events. 



o 

"5 0.8- 
o 

E 0.6- 



o 



0.4 -i 



0.2 f 



hi < 1.2 

- Data 

- Simulation 



ClulS, N/s = 7TeV 



<B 0.8- 

c 

o 

E 0.6 



3 4 567 10 20 30 

muon (GeV/c) 



0.4- 



0.2 s 



^B-tiWil-tSl3MH=| 



100 



hi < 1.2 

- Data 

- Sitnuiation 



ClulS, N/s = 7TeV 



o 

05 0.8- 

c 

o 

3 

E 0.6- 



3 4 567 10 20 30 

muon 'p^ (GeV/c) 



100 



0.4- 
0.2- 

oti^ 



hi < 1.2 

- Data 

- Simulation 



CMS, \/s = 7TeV 



3 4 567 10 20 30 100 

muon 'p^ (GeV/c) 



g 

a> 0.8 
o 

E 0.6 
o 

" 0.4 



0.2 



1.2 < hi <2-^ 
■ Data 
s- Simulation 



CMS, \/s = 7TeV 



2 3 45 



10 



. _ 20 30 

muon p^ (GeV/c) 



1- 




<U 0.8- 

§ - 

i 0.6=*^ 

LL 
CL 

0.4- 
0.2- 



100 



2 3 



1.2 < hi « 2.4 
• Data 

s- Simulation 



CMS, s/s = 7TeV 



10 



. _ 20 30 

muon p^ (GeV/c) 



1- 



o 
c 

03 
O 
3= 

<U 0.8- 

c 

o 

E 0.6- 



0.4- 



0.2- 



100 



1.2 < hi < 2.4 
• Data 
s- Simulation 



2 3 45 



CMS, s/s = 7TeV \ 

10 20 30 100 

muon p^ (GeV/c) 



Figure 11: Tag-and-probe results for the muon efficiency ^rec+id in data compared to simulation. 
Given that a tracker track exists, the plots show the efficiency as a function of muon px for Soft 
Muons (left), Particle-Flow Muons (middle), and Tight Muons (right) in the barrel and overlap 
regions (top), and in the endcaps (bottom). The measurement is made using ]/ip — > V-^V^ 
events for pj < 20 GeV/c and Z — > events for px > 20 GeV/c. For px < 3 GeV/c, to reduce 

the background, only tracks with MIP signature are considered. 



The tag-and-probe results in data and in simulation agree within the statistical uncertainties of 
the measurement almost everywhere. The only significant discrepancy is in the barrel around 
the turn-on of the efficiency curves, where the efficiency in data is systematically higher than 
in the simulation. This discrepancy arises from a small difference in the widths of the track-to- 
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Table 3: Muon efficiencies at the efficiency plateau for the different muon selections: efficiency 
measured from data, and ratio between the measurements in data and simulation. The first 
uncertainty quoted on the scale factor is the uncertainty on the efficiencies in data and simu- 
lation from the fitting procedure, which includes the statistical uncertainty; the second is from 
the additional systematic uncertainties described later in this section. 



Muon selection 

Region 


Eff. [%] Data/Sim. ratio 


Z — > 

Eff. [%] Data/Sim. ratio 


Soft 0.0 <\r]\ < 1.2 
1.2 < \r]\ < 2.4 


98.4+°| 1.010 ± 0.003 ± 0.010 
98.0toj 1.002 ± 0.007 ± 0.014 


99.2+°;^ 1.014 ±0.001 ±0.002 
99.9^^;^ 1.005 ±0.002 ±0.004 


Particle- 0.0 < \r]\ < 1.2 
Flow 1.2 < \r]\ < 2.4 


98.8;{^| 0.993 ± 0.003 ± 0.010 
98.4+°;^ 0.988 ± 0.007 ± 0.014 


99.7toi 0.999 ± 0.001 ± 0.002 
99.8^^:^ 0.999 ± 0.002 ± 0.004 


Tight 0.0 <\r]\ < 1.2 
1.2 < \r]\ < 2.4 


98.4+°J 0.998 ± 0.004 ± 0.010 
96.8+°;^ 0.979 ± 0.007 ± 0.014 


96.4+° 2 0.999 ± 0.002 ± 0.002 
96.0+° 3 0.983 ± 0.003 ± 0.004 



segment pulls in data and in simulation discussed in Section [43] the efficiency of the track-to- 
segment matching is slightly higher in data, and in the region of rapidly rising efficiency the 
effect is amplified by the large variation of the efficiency in the bin. The 1-2% data-simulation 
difference in efficiency for Tight Muons in the endcaps is explained by the fact that several 
CSCs not operational during most of the 2010 data taking were simulated as fully efficient; 
this has a negligible effect on efficiencies of the other muon selections because they require a 
match with only a single muon segment. Using a small sample of simulated events, we have 
verified that when these chambers are properly accounted for in the simulation, the efficiencies 
for Tight Muons in data and in simulation agree to better than 1%. 

The efficiency for the Tight Selection measured on muons from Z — > is slightly lower 

than that measured on muons from J/t/j — > ■ This difference is partly due to the bias 

introduced by the muon-plus-track trigger on the track quality criteria described previously 
and partly due to the different kinematics of the probes. The effect is well reproduced by the 
simulation. 

For Soft Muons and Particle-Flow Muons the plateau of the efficiency is reached at px ~4 GeV/c 
in the endcaps and ~6 GeV/c in the barrel, while for Tight Muons it is reached at ~10 GeV/c in 
both regions. The values of efficiencies at the plateau region obtained using ]/ip — )■ and 
Z — )■ events in data and simulation are given in Table |3j The efficiencies are high, and 

data and simulation are generally in good agreement. The plateau efficiency for Soft Muons 
is 1-2% higher in data than in simulation, again due to a difference in the widths of the track- 
to-segment pulls. The Particle-Flow and Tight Muon selections are much less affected by this 
difference because they use looser matching criteria between tracks and muon segments. The 
efficiency at the plateau for the Particle-Flow Muon selection is very close to 100% because the 
algorithm applies relaxed selection criteria to the high-pj muon candidates if they are isolated. 

The dependency of the plateau efficiency on the pseudorapidity is measured using Z — > f/+f<^ 



events and is shown in Fig. 12 The data and simulation agree to better than 2%. 



To estimate the effect of pile-up on muon identification performance, the efficiency at the 
plateau is measured for both ]/\p — > ji^ji^ and Z — > events as a function of the num- 

ber of reconstructed primary vertices. No loss of efficiency is observed for events containing 
up to six reconstructed primary vertices, the maximum multiplicity for which a measurement 
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Figure 12: Muon efficiency erec+id in data and simulation as a function of muon pseudorapidity 
for Soft Muons (left), Particle-Flow Muons (middle), and Tight Muons (right). The efficiencies 
were calculated relative to the tracker tracks with pj > 20 GeV/ c by applying the tag-and-probe 
technique to Z — > events. 



could be made with a statistical uncertainty below 10%. 

Trigger efficiencies obtained by a similar tag-and-probe technique are described in Section |9] 
5.1 .3 Systematic uncertainties 

The various contributions to the possible systematic bias in the data-to-simulation ratios of 
efficiencies calculated using the tag-and-probe method are estimated using simulated and real 
data. 

Bias in the measured efficiencies that could be introduced by the tag-and-probe method and 
its implementation is studied by comparing the efficiencies obtained by applying the tag-and- 
probe method to simulated data containing — >• (Z — )> //+f/^) decays and various back- 
ground contributions with the "true" efficiencies computed by simple counting of the passing 
and failing probes in "j/xp — > (Z — ?► ]d^ji^) MC events. The difference in the efficiencies 
is less than 0.5% for muons from J/t/> — ?► It is also less than 0.5% for most muons from 
Z fi^fi' (20 < Pt < 60 GeV/c), and less than 1.5% for the others. The differences are compat- 
ible with zero within the statistical uncertainties; hence, no systematic uncertainty is assigned 
as a result of this test. 

For ]/xp events, the efficiencies are recomputed with the tag-and-probe method using a simple 
Gaussian instead of a Crystal Ball function [28J to model the resonance and with a quadratic 
polynomial instead of an exponential to model the background. The differences in the efficien- 
cies resulting from this variation in the assumed signal shape are under 0.1%. The efficiencies 
obtained with a polynomial background shape are systematically ^1% higher than those ob- 
tained using an exponential background. The difference between the two results is taken as a 
conservative estimate of systematic uncertainty in the background modelling. The same effi- 
ciencies have also been recomputed without the requirement that the probe tracks have a MIP 
signature in the calorimeters. The results are fully compatible with those obtained with a MIP 
requirement, but have larger uncertainties. Simulation shows that in this low-pj range tag- 
and-probe efficiencies estimated with a MIP requirement are systematically higher, by 1-2%, 
than without a MIP requirement, due to small correlations between the energy deposition in 
the calorimeters and the number of hits in the muon chambers. This bias cancels out in the 
data/ simulation ratio, so no corrections accounting for it are made. 
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For Z —> events, the efficiencies are recomputed using only isolated probe tracks, which 

reduces the background by a factor of two. The results agree with those from all probes at the 
level of 0.1%. As a conservative estimate of systematic uncertainty on the plateau efficiencies 
resulting from the background estimation, the scale of the largest difference between this es- 
timate and that from simulation is taken; this amounts to 0.2% in the barrel and 0.4% in the 
endcaps. 

For ]/xp events, the kinematic distributions of the signal probes are extracted from the distri- 
butions of all probes by using the SPlot technique [29|. The distributions were found to be in 
good agreement with those predicted by simulation, and therefore no systematic uncertainty is 
assigned to the procedure of averaging the efficiencies for different probes within each {pt,i]) 
bin. 

Possible bias in the measurements of single-muon efficiencies due to the presence of a second, 
tag muon in the event is studied by changing the separation criteria from the angular sepa- 
ration AR > 0.5 to the tighter requirement that the distance between the coordinates of the 
two muons in the innermost muon station be larger than 2 m and by using only the pairs of 
muons that bend away from each other in the magnetic field inside the solenoid. The effect 
on the efficiencies measured using ]/xp events is 1% in the endcaps and 0.3% in the barrel; this 
difference is taken as an additional systematic uncertainty on the data-to-simulation ratios of 
efficiencies. The impact on the efficiencies measured using Z — > events is negligible, and 

no additional systematic uncertainty is assigned to them. 

In the measurements of efficiencies using ]/\p — > events, no attempt is made to separate 

promptly produced J/t/J from those originating from the decay of b quarks. Differences in 
efficiencies obtained using these two samples of muons are studied in simulation and have 
been found to be less than 1% and compatible with the statistical uncertainties on each. In the 
kinematic range over which the measurement is made, the fraction of J/i/?'s from the decays of 
b quarks is always below 50% [18]/ so the possible effects on the data-to-simulation efficiency 
ratios are below 0.5%. No additional systematic uncertainty is assigned. 

The various contributions to systematic uncertainty were combined in quadrature; the overall 
uncertainties at the efficiency plateau are shown in Table |3j 

5.2 Reconstruction and identification efficiency for nearby muons 

If a muon has one or more other muons in its vicinity, their signals in the muon system could 
overlap, resulting in identification efficiencies lower than for single or well-separated muons. 
For example, such topologies are common for muon pairs produced in the decays of low-mass 
resonances such as ]/'ip. Another example is hypothetical highly collimated leptons, also re- 
ferred to as "lepton jets", predicted in different models ||30H32| proposed to explain the excess 
of cosmic-ray leptons in recent astrophysical observations P33I - I35I . In this section, we report 
the measurement of reconstruction and identification efficiency for such nearby muons. 

The muon identification performance for nearby muons is studied using the data collected 
during the 2010 LHC data-taking period and compared to the expectations from simulation, 
in which boosted muon pairs are generated using PYTHIA (S). Two muon selections are con- 
sidered: Tracker Muons with at least two tightly matched segments in the muon system, and 
Tight Muons. Studies on simulated events have shown that the purity of Tracker Muons with 
two or more well-matched muon segments is similar to that of Tight Muons. 

In this study, the efficiency of identifying nearby muons as Tight or Tracker Muons is mea- 
sured using a sample of dimuons from the decays of low-mass resonances: ]/ip, (p, and p/ co. 



5.2 Reconstruction and identification efficiency for nearby muons 
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These resonances provide muon pairs with kinematic and topological properties similar to 
those expected for hypothetical collimated muons, notably a small angular separation between 
the muons. 

The sample used for this study consists of pairs of tracker tracks each with pj above 5 GeV/ c, 
associated with the same primary vertex, and for which the invariant mass is in the vicinity 
of the invariant mass of one of the above resonances. Because the lower mass resonances {(p 
and p/co) have very large combinatorial background, both tracks are also required to have one 
loosely matched segment in the muon system. To check for any bias introduced by this re- 
quirement a measurement was made with tracker tracks from J/i/? decays identified as possible 
muon candidates by using only the calorimeter information. This gives very similar results 
hence showing any bias is small. The contribution from the residual background is evaluated 
from fits to the side bands, by a procedure similar to that used in the tag-and-probe method 
(Section |5.1| , and suitably subtracted. The efficiency is defined as the ratio of the number of 
dimuon candidates passing both the above selection and the muon identification under study 
to the number of dimuons passing the above selection. 
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Figure 13: Efficiency for identifying both muons in the dimuon pair as Tight (left) and Tracker 
(right) Muons as a function of the angular separation of the two tracks computed at the surface 
of the first muon station. Measurements obtained using ]/tp (squares), (p (inverted triangles), 
and p/co (triangles) are compared with the expectations from the simulation (circles). 



The efficiencies for identifying both muons in the dimuon pair as Tight or Tracker Muons are 
shown in Fig. [l3] The efficiencies are plotted as a function of the angular separation of the 
two tracks, AR = ^y (Arj)^ + (Acp)^, computed at the surface of the first muon station; the 
collimated muons are expected to populate the range of AR < 0.7. For Tight Muons, a drop 
in efficiency at small values of AR is observed; this inefficiency is introduced by a cleaning 
procedure used at the seeding stage of the global muon reconstruction to eliminate muon seeds 
leading to duplicate muons. The efficiency for Tracker Muons, however, remains high at all 
AR values, demonstrating that Tracker Muons are fully adequate for studies involving nearby 
muons. For both types of selections, the results obtained using different resonances are in good 
agreement, demonstrating that the dependence of the efficiency on AR is not affected by the 
decay kinematics and combinatorial background. The results of the measurements are also 
well reproduced by the Monte Carlo simulation. 
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5.3 Muon identification probability for particles other than muons 

One can obtain pure samples of kaons, pions, and protons from resonances of particle decays 
such as Kg — > TT+TT^, A — > pn^ (and charge conjugate), and (p — ?> K+K^. The resonances are 
reconstructed using pairs of tracker tracks that are associated with a common decay vertex, 
with a selection similar to that described in Ref. [36 J. In A decays, the highest momentum track 
is assumed to be that of the proton. A data sample collected with a jet trigger (minimum pj of 
15GeV/c) is used, and simulated QCD events, filtered using the same jet trigger, are used for 
comparison. The simulated events have been reweighted to account for a small difference in 
the hadron momentum spectrum with respect to the data sample. 

We compute the fraction of events in which a hadron track is identified as a Soft Muon, Particle- 
Flow Muon, or Tight Muon as a function of several relevant track parameters. Background 
subtraction using resonance sidebands is performed to determine the muon misidentification 
probability for the particles under study. Invariant mass spectra are fit with a sum of signal 
and background shapes, using a double Gaussian for the signal and a power law for the back- 
ground. One fit to the entire mass spectrum is made for each resonance to provide the scale 
factor between the number of hadrons counted in the sideband region and the background 
estimation in the signal region. The scaled number of hadrons in the sideband region is then 
subtracted from the number counted in the signal region, in each bin of the distribution of 
the hadron track parameter under study. The same background-subtraction procedure is then 
repeated only for hadrons that share the tracker track with that of a muon. By dividing the 
sideband-subtracted number of muon-matched hadrons by the sideband-subtracted number 
of hadrons before any matching to muons, we obtain a misidentification probability for a given 
hadron type. The same method is applied to events in data and simulation. 

In addition to punch-through and decay in flight, there is a third mechanism by which hadrons 
can be misidentified as prompt muons. This mechanism is random matching between the 
hadron track in the inner tracker and a track stub in the muon system from one of the other 
tracks in the jet that may be due to a muon. The frequency of random matching is sensitive to 
the particular event topology of the sample. For example, in jet-triggered events, the increase in 
the average number of tracks per event in comparison to minimum-bias-triggered events leads 
to an increased probability of random matching. To illustrate the effect of random matching. 



we present the proton-to-muon misidentification probabilities as a function of Njracks in Fig- 14 
where Njracks is the number of tracks in the vicinity of the proton track, within a cone of radius 
AR < 0.2. It is clear for both data and simulation that the misidentification probability increases 
with Njracks especially in the Soft Muon selection. To remove much of the contribution due to 
random matching from the probability of misidentifying a hadron as a muon, we impose a 
requirement of Njracks < 4 for the rest of the results in this section. 



The resulting muon misidentification probabilities are shown in Figs. 15 and 16 as a function 
of particle momentum and pseudorapidity respectively. The shapes of the distributions, well 
reproduced by simulation, are due to a combination of acceptance (a minimum momentum is 
required to reach the muon system), the amount of material before the muon system, and the 
distance available for pions and kaons to decay before reaching the calorimeters. For pions and 
kaons, the misidentification probabilities are below 1% for all muon selections and decrease at 
p > 10-15 GeV/c due to fewer of the hadrons decaying to muons within the detector volume. 
For protons, the probability to be identified as a muon slowly increases with momentum but 
remains low in the accessible momentum range, which confirms that punch-through is small 
and that at low momenta the main reason for misidentification of pions and kaons is decays in 
flight, in agreement with the predictions from simulation discussed in Section |4j As expected. 
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Figure 14: The fraction of protons that are misidentified as a Soft Muon (left), Particle-Flow 
Muon (centre), or Tight Muon (right) as a function of Nxracks/ where Njracks represents the num- 
ber of tracks in the vicinity of the proton track (with AR = ^J {Ar])^ + {A(p)^ < 0.2). Only 
protons with p > 3GeV/c are included. The first bin includes events with Nxracks = srid 1, 
the second bin includes those with Njracks = 2 and 3, etc. The uncertainties indicated by the 
error bars (data) and shaded boxes (pythia simulation) are statistical only. Negative values 
arise from statistical fluctuations in the number of events in the signal and sideband regions. 



the misidentification probabilities are found to be independent of the azimuthal angle and the 
decay length of the mother particle within the statistical uncertainty. 

Overall, the probability to misidentify a hadron as a muon is the largest for Soft Muons, de- 
creases slightly for Particle-Flow Muons, and drops significantly for Tight Muons. As shown in 



Section 5.1 the lower misidentification probability for Tight Muon selection comes at the cost of 
a few percent lower muon identification efficiency. It is this trade-off between misidentification 
probability and efficiency that motivates using different muon selections for different analyses. 



6 Muon Momentum Scale and Resolution 

The measurement of the muon transverse momentum is highly sensitive to the alignment of the 
tracker and of the muon chambers, to the composition of material and its distribution inside the 
tracking volume, and to the knowledge of the magnetic field inside and outside the solenoid 
volume. 

The relative bias A(px)/pT in reconstructed muon transverse momentum with respect to its 
true value that could be caused by imperfect knowledge of the magnetic field is generally con- 
stant as a function of momentum. Similarly, inaccuracies in the modelling of the energy loss 
(dependent on the material distribution) produce relative biases that are essentially indepen- 
dent of the muon momentum. On the other hand, alignment effects produce relative biases 
that generally increase linearly with momentum. 

The momentum scale and resolution of muons are studied using different approaches in dif- 
ferent Pt ranges. At low and intermediate px (^lOOGeV/c), the mass constraint of dimuon 
decays from the ]/xp and Z resonances is used to calibrate the momentum scale and measure 
the momentum resolution. In the high-px range (^100 GeV/c), the muon momentum scale and 
resolution can be measured using cosmic-ray muons (with the exception of the high- 1 tj \ region) . 

The lower px range of the muon spectrum, px < 10 GeV/c, has been studied in Ref. [37]. In this 
region of px, alignment effects are less important, and biases in the reconstructed momentum 
mostly arise from uncertainty in the modelling of the detector material and in the description of 
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Figure 15: The fractions of pions (top), kaons (centie), and protons (bottom) that are misidenti- 
fied as Soft Muons (left), Particle-Flow Muons (centre), or Tight Muons (right) as a function of 
momentum. Only particles with Njracks < 4 are included. The uncertainties indicated by the 
error bars (data) and shaded boxes (pythia simulation) are statistical only. 



the magnetic field used when reconstructing the track. Results obtained using ]/xp events show 
that the overall relative bias in the tracker measurement of the muon pj in this momentum 
range is ^0.1%. The muon pj resolution cr{pj)/pj was found to be between 0.8% and 3% 
depending on and in good agreement with the simulation. 

In the intermediate-px range, two approaches to study the muon pj measurement have been 
developed. The first, referred to as MuScleFit (Muon momentum Scale calibration Fit), pro- 
duces an absolute measurement of momentum scale and resolution by using a reference model 
of the generated Z lineshape convoluted with a Gaussian function. The second, called SIDRA 
(Simulation DRiven Analysis), compares the data with the full simulation of the Z decay to 
two muons in the CMS detector and provides a way to directly modify the simulation to better 
match the data. As these two methods have different approaches to the same problem, the 
difference between the results provides a useful crosscheck and gives an estimate of a system- 
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Figure 16: The fractions of pions (top), kaons (centre), and protons (bottom) that are misiden- 
tified as Soft Muons (left), Particle-Flow Muons (centre), or Tight Muons (right) as a function 
of pseudorapidity. Only particles with p > 3GeV/c and Njracks < 4 are included. The uncer- 
tainties indicated by the error bars (data) and shaded boxes (pythia simulation) are statistical 
only. 



atic uncertainty in the measurement. The results obtained with these methods are reported in 
Section l6m 

At high Pj (Section [6.2| , the resolution is determined by comparing cosmic-muon tracks recon- 
structed independently in the upper and lower halves of the detector, while the scale bias is 
evaluated by using what is called the "cosmics endpoint method". 

6.1 Measurements at intermediate pj 

As previously mentioned, a sample of muons produced in the decays of Z bosons is well suited 
for measuring the muon momentum scale and resolution in the intermediate range of trans- 
verse momentum, 20 < px < lOOGeV/c. Muons from Z-boson decays are identified using the 
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Figure 17: The position of the Z peak reconstructed by the MuScleFit method as a function of 
muon (p for positively charged muons (left), (p for negatively charged muons (middle), and r] 
for muons of both charges (right). Results obtained from data (black triangles) are compared 
with those from simulation (red circles). The (p plots also show the results of sinusoidal fits; the 
values of the fit parameters are given in the text boxes below the labels. 

Tight Muon selection. They are compared to the simulated sample of Z — )■ and Drell- 

Yan dimuon events reconstructed using the realistic misalignment scenario (see Section |2]) and 
passing the same set of muon selection criteria. Tight Muon selection is chosen since it is used 
in all electro-weak precision measurements in CMS. However using a different selection has 
no significant impact on the results presented in this section. 



Measurements of the muon momentum scale and resolution, reported in Sections 6.1.1 and 



6.1.2 refer to the px of the tracker-only fit as it is the default muon momentum assignment in 
the range of px < 200 GeV/c (see Section|3]|. Studies of the resolution of the px measured by the 
muon system only, provided by the standalone-muon track fit, are described in Section 6.1.3 



6.1.1 Muon momentum scale 

In the MuScleFit approach, the biases in the reconstructed muon px are determined from the 



position of the Z mass peak as a function of muon kinematic variables. Figure 17 shows the 
position of the peak determined from a fit to a Voigtian (the convolution of a Lorentzian and a 
Gaussian) in bins of muon (p, shown separately for positively and negatively charged muons, 
and muon rj, for both data and simulatioij^ We observe a sinusoidal bias as a function of 
(p, antisymmetric for muons of opposite charges. The maximum shift of the position of the 
mass peak is about 0.5%. This is significantly larger than the relative bias observed at low 
Px using the ]/\p events, indicating that the relative bias is px dependent and that its main 
source at intermediate px is the residual tracker misalignment. As explained in Section |2| the 
misalignment scenario used in the simulation was obtained by the alignment procedure similar 
to that used in the data but was based only on the sample of cosmic-ray muons available at the 
time of the MC production. This results in a different phase and slightly larger amplitude of the 
bias in MC simulation. In addition, a small /^-dependent bias is also present. This bias exhibits 
no dependence on muon charge, has a parabolic shape, and is reasonably well reproduced by 
the Monte Carlo simulation. 

MuScleFit uses an unbinned likelihood fit with a reference model to correct the momentum 
scale. Given the shape of the Z-peak position biases, the following ansatz function is used for 
the calibration of the momentum scale: 



^Note that the peak position returned by the fit is not expected to perfectly match the PDG value of the Z mass. 
The generator-level Z lineshape is not symmetric around the peak and has a higher tail at lower mass values. When 
it is convoluted with the detector resolution effects (which can be approximated with a Gaussian), the peak shifts 
towards lower values. A simplified test performed by convolving the reference lineshape with a Gaussian with 
a = 1.5 GeV/c^ gives a peak position of 90.8 GeV/c^. 
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Figure 18: The position of the Z peak in data as a function of muon cp for positively charged 
muons (left), (p for negatively charged muons (middle), and r] for muons of both charges (right) 
before (black circles) and after (red triangles) the MuScleFit calibration. The (p plots also show 
the results of sinusoidal fits; the values of the fit parameters are given in the text boxes below 
the labels. 



Pj = Pt(1 + b ■ Pj + c ■ t] + q ■ d ■ Pj ■ sin(^ + e)) , (2) 

where q is +1 for fi^ and —1 for fi^ , and b, c, d, and e are the fit parameters. The model taken as 
calibration reference is a lineshape of the Z decaying to dimuon pairs as described in Ref. 
and generated with a high granularity in 1001 bins between 71.2 and 111.2 GeV/c-^. Figure 
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shows the results of this calibration procedure on data. The calibrated position of the mass peak 
is consistent with being flat within the statistical uncertainties, demonstrating that the biases 
are successfully removed. When averaged over <p, the correction is small, and the position of the 
dimuon invariant-mass peak remains practically unchanged. The same calibration procedure 
also successfully eliminates momentum scale biases present in the simulation. 

The strategy implemented in the SIDRA method consists of modifying the reconstructed inva- 
riant-mass spectrum of simulated Z — )■ events by additional shifts and resolution distor- 
tions to make it agree with the invariant-mass distribution observed in data. This means that 
unlike the MuScleFit method, the SIDRA method calibrates only relative biases between data 
and simulation. This approach assumes that the resolution in data is slightly worse than in 
simulation and is well suited for the present study. The Z — )■ fi^ji^ candidates in data and in 
simulation are binned according to their reconstructed parameters, and the difference between 
the two distributions is minimized by using a binned maximum likelihood fit. At each min- 
imization step, the reconstructed transverse momentum pj of the simulated muons, px, sim/ is 
modified as follows: 



— = ^ Sf:^{q,(p,r]) + aKj{q,(p,}]) Gauss{0,l) , (3) 

where S,(j and (Tkj are parameters controlling the scale shifts and resolution distortions, re- 
spectively, and Gauss (0, 1) denotes a sampling according to a Gaussian of zero mean and unit 
variance. The fit parameters depend in general on the muon charge q, its azimuthal angle (p, 
and its pseudorapidity rj. This ansatz assumes that the differences between the data and the 
simulation are due to misalignment, with relative effects increasing with pj. This assumption 
is justified by the excellent agreement between data and simulation for low-mass dimuon reso- 
nances [i37il . Similarly to the MuScleFit case, several exploratory studies suggest dependencies 
of the type 

SKj{q,(p,f]) = A + Bt]^ + q C sm{(p - (po) ; (4) 
crKj{q,(p,v) = A' + B'f]^, (5) 
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where A, B, C, (po, A', and B' are the parameters to be determined in the fit. The dependence on 
charge iniluences the choice of the binning for the fit: since the proposed ansatz function has a 
charge-dependent term as a function of (p, we employ a two-dimensional grid, binning events 
according to the reconstructed dimuon invariant mass and the azimuthal angle of one of the 
two muons. 
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Figure 19: Top: distributions of the dimuon invariant mass for the selected Z — > can- 
didates in data (points with error bars) and in simulation without ("reference MC") and with 
("corrected MC") corrections from SIDRA applied. Bottom: bin-by-bin difference (rebinned for 
clarity) between the simulation and the data, divided by the expected statistical uncertainty, 
for MC samples without (filled black circles) and with (open red circles) the SIDRA corrections. 
The uncertainties are statistical only. 
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The results of the application of the SIDRA method to the dimuon invariant-mass spectrum 



in 2010 data are shown in Fig. 19 Applying SIDRA corrections to the simulation improves 



agreement with the data. The scale shifts Skj as a function of cp and rj are shown in Fig. 20 



superimposed with the corresponding shifts obtained with the MuScleFit method. As also 



shown in Fig. 17 the phases and amplitudes of the (^-dependent biases present in data and 
in simulation are different: for px ~ / 2, characteristic of this study, the amplitude of the 
sinusoidal correction to be applied to the simulation in order to obtain the best match with the 
data is ?»1.5%; differences between the corrections from the two methods do not exceed 0.3%. 
When examined as a function of tj, the scale shifts between data and simulation are consistent 
with zero in the barrel region and increase with \r]\. For pi = Mz/2 and \rj\ > 2, they are 
?»0.5% with a difference between the two methods of ^0.1%. When integrated over (p and tj, 
the overall difference in muon momentum scale Skj between data and simulation is found to 
be 0.016 ± 0.012 (stat.) c/TeV for SIDRA and 0.020 ± 0.006 (stat.) c/TeV for MuScleFit. 

We use the difference between the MuScleFit and SIDRA relative simulation-to-data corrections 
as the systematic uncertainty in the measurements of the absolute momentum scale bias ob- 
tained with MuScleFit. Another source of systematic uncertainties are the theoretical uncertain- 
ties in the reference models used by both SIDRA and MuScleFit, such as uncertainties in parton 
distribution functions, initial state radiation and higher-order QCD effects, and weak and QED 
interference effects. Dedicated simulation studies show that these uncertainties can produce 
shifts in the Z mass peak position of at most 0.1%, which corresponds to an uncertainty in pj 
scale of order 0.1% ■ px/ [91 GeV/c] if positive and negative muons are equally affected. Sum- 
ming both types of systematic uncertainties in quadrature, the amplitude of the (/^-dependent 
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Figure 20: Comparison of the differences between muon momentum scale in data and that in 
the simulation obtained with the SIDRA and MuScleFit methods as a function of the azimuthal 
angle for positive and negative muons (left) and pseudorapidity (right). Only statistical uncer- 
tainties are shown. 



scale correction at px = Mz/2 is found to be 0.266 ± 0.010 (stat.) ± 0.046 (syst.)GeV/c. When 
integrating over (p, the scale correction for the same pj varies from A(px) = —0.130 ± 0.022 
(stat.) ± 0.046 (syst.) GeV/c for f/ = to A(px) = 0.234 ± 0.048 (stat.) ± 0.046 (syst.) GeV/c for 
= 2.1. The future versions of the alignment workflow will include a Z-mass constraint as an 
integral part of the alignment procedure, which will strongly reduce the (p and t] dependence 
of the momentum bias. Averaged over the whole acceptance, the relative bias in the muon 
momentum scale is measured with a precision of better than 0.2% and is found to be consistent 
with zero up to pj values of 100 GeV/c. 

6.1.2 Muon momentum resolution 

The techniques used to calibrate the muon momentum scale can also be used to measure the 
muon transverse momentum resolution. In SIDRA, an extra smearing in data with respect to 
the resolution predicted by the simulation is obtained from the parameter, as explained in 



Section 6.1.1 The full resolution in data can then be evaluated by applying the scale shifts and 
extra smearing to the reconstructed pj of muons in simulation and comparing the obtained 
Pj values with the "true" pj. The MuScleFit procedure takes into account the correlations 
between the two muons, but neglects contributions to the mass resolution from the (p and ?/ 
resolutions, as well as covariance terms. These contributions were found to be small with 
respect to other terms at the typical pj range of muons from Z; nevertheless, the systematic 
uncertainty due to these approximations is estimated from the simulation and accounted for in 
the final uncertainty of the result. Resolution and scale biases are fit in successive iterations to 
minimize correlation effects. The assumed functional form to fit the momentum resolution in 
data contains a term linear in pj, and separate terms for positive and negative t], each parabolic 
in T] with a common minimum at ?/ = 0. In the simulation many more events are available, and 
the 7] component contains a third term, symmetric about rj = 0. 



The calibration of the momentum scale described in Section 6.1.1 improves the muon momen 



tum resolution by ^2%. Figure 21 shows the MuScleFit and SIDRA measurements of muon 
transverse momentum resolution versus ?/ after correcting for biases in the momentum scale, 
for both data and simulation. The grey band enveloping MuScleFit results in data shows sta- 
tistical and systematic uncertainties summed in quadrature. The uncertainties from the choice 
of a particular function for the resolution shape and the approximations in the method are es- 
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Figure 21: Relative transverse momentum resolution cr['pj)/pj in data and simulation mea- 
sured by applying the MuScleFit and SIDRA methods to muons produced in the decays of Z 
bosons and passing the Tight Muon selection. The thin line shows the result of MuScleFit on 
data, with the grey band representing the overall (statistical and systematic) la uncertainty 
of the measurement. The circles are the result of MuScleFit on simulation. The downward- 
pointing and upward-pointing triangles are the results from SIDRA obtained on data and simu- 
lation, respectively; the resolution in simulation was evaluated by comparing the reconstructed 
and "true" px once the reconstructed pj was corrected for (^-dependent biases. The uncertain- 
ties for SIDRA are statistical only and are smaller than the marker size. 



timated by comparing the result of the fit using the same function in simulation with the true 
MC resolution. The bin-by-bin difference between the two results is taken as the systematic 
uncertainty. The relative difference is on average 6% with an RMS of 4%. Another source of 
systematic uncertainty included in the band is theoretical uncertainties in the reference models 



discussed in Section 6.1.1[ They can produce an extra smearing of the Z mass distribution of at 
most 0.5%, which can be interpreted as an uncertainty of 0.5% / \/2 for muons with pj Mz / 2. 
The statistical and systematic uncertainties are of similar magnitude; the overall (statistical and 
systematic combined in quadrature) la uncertainty of the measurement varies from 20% to 
40% of the resolution in the studied acceptance range. 



As can be seen in Fig. 21 the results obtained with the two methods agree within the uncer- 
tainties: the largest difference in the barrel is {a{pj)/ pT)MuScieFit — (c^(Pt)/Pt)sidra = 0.003 ± 
0.003(stat. syst.), while in the endcaps it is —0.018 ± 0.013(stat. © syst.). The relative pj res- 
olution in the intermediate px range obtained using MuScleFit is found to be in the range from 
1.3% to 2.0% for muons in the barrel and up to ^6% for muons in the endcaps, in good agree- 
ment with the results obtained from simulation. The a{pj) / pj averaged over (p and rj varies in 
Px from (1.8 ± 0.3(stat.))% at px = 30GeV/c to (2.3 ± 0.3(stat.))% at px = 50GeV/c, again in 
good agreement with the expectations from simulation. 

6.1 .3 Momentum resolution of standalone muons 

The momentum resolution for standalone muons is estimated using the px of the tracker track 
as reference: 
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Rsta(l/PT)= ((l/pT)^'^-(l/pTrM /(1/Pt) 



\trk 



(6) 



As the resolution of the tracker tracks at low and intermediate pj is expected to be about an 
order of magnitude better than the resolution of the standalone-muon tracks, Eq. (|6| provides 
a good estimate of the resolution for standalone muons. The relative difference between the 
resolutions measured with respect to the tracker-track pj and the true pj was evaluated from 
simulation and found to be smaller than 1% in the barrel and smaller than 5% in the endcaps. 

The sample of muons used for this study is selected from events passing the single-muon trig- 
ger with a minimum pj threshold of 5 GeV/ c by applying the standard requirements for Tight 
Muons (see Section |3]). Selected muons were subdivided into subsamples according to the pj 
and ?] of the tracker track, and a Gaussian fit to the distribution of Rgta (Eq. (|6]l) is performed for 
each subsample, with the fit restricted to a range ±1 RMS about the sample mean. The resolu- 
tion thus measured using data is compared with results obtained from simulation. Simulated 
muons are selected from MC samples of QCD, W — > j-iv, and Z — > ji^ji^ events using the same 
trigger and selection criteria as applied to data. 
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Figure 22: The l/px resolution of standalone-muon tracks with respect to tracker tracks, as a 
function of rj of the tracker track for muons with pj > 15 GeV/c. Resolution is estimated as the 
(7" of a Gaussian fit to Rgta (1 / pj) defined in Eq. (|6j|. 



We report results for standalone muons with a beam-spot constraint applied in the track fit. 



Figure 22 shows the widths of the Gaussian fits to the distributions of Rgta as a function of rj for 
muons with pj > 15GeV/c. The resolution in the barrel remains better than 10% up to px = 
100 GeV/c. The overall shape of the tj dependence is reproduced by the MC simulation to within 
10-15%. Considering that standalone muons are usually not used in physics analyses on their 
own and mainly serve as a component of the global-muon reconstruction, the performance of 
which is well described by the simulation, this difference between resolutions for standalone 
muons in data and simulation is acceptable. A slight asymmetry between resolutions in the 
negative and positive endcaps is due to small differences in the precision of the alignment of 
the muon chambers. The bias in the momentum scale of standalone muons, given by the mean 
values of the fits to the Rgta distributions, does not exceed 1% in the barrel region and 5% in the 
region of 0.9 < |?/| < 2.1. 
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6.2 Measurements at high 

High-px muons are an important signature in many searches for new physics, so it is crucial 
that the performance of their reconstruction, which has some significant differences to that of 
lower-px muons (such as an increased role of the muon system in momentum measurement 
and larger impact of showering) is well understood. While there are few high-px muons from 
LHC collisions in the 2010 dataset, cosmic-ray events provide a source of muons with a momen- 
tum spectrum extending to quite high px; cosmic-ray muons having momentum up to a few 
TeV/ c have been observed in CMS. Here we present results using cosmic-ray events collected 
during periods in 2010 when the LHC was not delivering collisions. 

6.2.1 Muon momentum resolution from cosmic-ray muon data 

Cosmic-ray muons that traverse the entire CMS detector can be used to evaluate the momen- 
tum resolution by comparing the momenta reconstructed independently in the upper and 
lower halves of the muon system, a procedure that was first applied to cosmic-ray muons 
collected in 2008 [2J. The angular distribution of cosmic-ray muons that traverse the detector is 
strongly peaked in the vertical direction. To select events in which the muons are most similar 
to those produced in collision events and used in physics analyses, we only use pairs of tracks 
that pass close to the interaction point, enforced in practice by requiring each track of the pair 
to contain hits from at least eight different layers of the silicon strip tracker and at least one 
hit from the pixel detector. The results presented in this section are for the barrel region only, 
because few cosmic muons cross both endcaps and pass near the centre of CMS. 

To estimate the muon q/pj resolution, we define the relative residual. 



(^/p ^)"PPe^_(^ /p^)lower 



where "upper" and "lower" refer to the results of the fits in the two halves of CMS, and the fac- 
tor of \/2 accounts for the fact that the two fits are independent. We also define the normalized 
residual (pull) as 



to examine the behavior of the uncertainties from the track fit. To study the various effects 
as a function of px, we examine both the truncated sample RMS and a Gaussian fit to the 
distributions of R(ij/ px) and P{ci/ pj), with the fit range in each px bin under consideration 
restricted to the region of ±1.5 RMS centred on the sample mean. The px used in the binning 
is that of the lower tracker-only track. 



Figure 23 'a) shows the Gaussian widths of the pulls as a function of px for the tracker-only and 
global fits, and for the sigma-switch and Tune P algorithms (see Section |3]|. While the cores 
of the pull distributions are well described by Gaussians, the tails are non-Gaussian, so these 
widths are slightly smaller than the sample RMS values, which are not shown here. The widths 
of the pulls for the tracker-only fit are within ±10% of unity over the entire px range studied, 
indicating that the track uncertainties are well estimated. The widths of the pulls for the global 
fit are close to unity in the range of px ^ 200 GeV/c, where the tracker information dominates 
in the fit, and start to become larger at px ~ 200 GeV/ c. One important factor contributing to 
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this is that the muon alignment position uncertainties, which account for the precision with 
which the positions of different detector components are known, are set to zero in the track 
fits. As a consequence, the total uncertainties are underestimated at high pj. As Tune P uses 
less information from the muon system than the global fit, the pulls are in between those of 
tracker-only and global fits. 

Figure |23|b) shows the sample means for the distributions of P{q/pj) for the same four muon 
reconstruction algorithms. The means are consistent with zero up to pj values of lOOGeV/c, 
above which fits with muon information begin to show a sizeable bias. The biases seen in the 
means of the residual distributions R{cj/ pj), which are not shown here, are small, indicating 
that the source of the bias in the means of P{q/pj) is most likely to be the underestimated 
uncertainties in the fit. This bias should become smaller as the muon system alignment and its 
uncertainties are further refined. 
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Figure 23: (a) Widths of Gaussian fits to the distributions of the muon q/pj pulls (as defined 
by Eq. (jsjl) for the tracker-only and global fits, and for the output of the sigma-switch and Tune 
P algorithms, as a function of the pj of the muon; (b) the means of the same distributions. 



Figures [24|a) and (b) show the relative resolution as measured by the Gaussian width and the 
truncated sample RMS, respectively, for the tracker-only and global fits, and for the sigma- 
switch and Tune P algorithms as a function of the pj of the muon. Table |4] summarizes the 
performance of these algorithms for muons with transverse momentum in the range 350 < 
Pt < 2000GeV/c. The Gaussian width gives a measure of the core resolution, while the trun- 
cated sample RMS includes the effects of the tails of the distribution; both can be separately 
important for considerations of momentum resolution, possibly depending on the details of 
the physics analysis being considered. 

Muon momentum resolution in the pj region below approximately 200 GeV/ c is dominated by 
the track measurement in the inner tracker, so one already has the best performance by using 
the tracker-only fit for low-to-intermediate-px rnuons. At high pj, the extended path length 
through the magnetic field between the tracker and the muon system leads to improved reso- 
lution, provided that the pattern recognition carefully selects hits in the muon system, avoiding 
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Figure 24: (a) Widths of Gaussian fits of the distributions of the muon q/pj relative residuals (as 
defined by Eq. ^) for the tracker-only and global fits, and for the output of the sigma-switch 
and Tune P algorithms, as a function of the pj of the muon; (b) sample RMS values (truncated 
at ±1) of the same distributions. 



Table 4: Comparison of the fitted Gaussian width a, sample RMS (truncated at ±1), and counts 
in the tails of the R{cj/ pj) distribution for 294 muons with measured transverse momentum in 
the range 350 < pi < 2000GeV/c. 



Fit/ selector 


Fitted a, % 


RMS, % 


R{q/PT) < -1 


R{q/pT) > 1 


Tracker-only 


7.3 ± 0.3 


9.7 ±0.4 








Global 


7.8 ± 0.3 


15.6 ± 0.7 


4 


9 


Sigma switch 


7.1 ± 0.3 


9.8 ± 0.4 








Tune P 


6.2 ± 0.3 


8.2 ± 0.3 









hits due to showers and hits in any muon chamber whose position is poorly known. In Fig. 24 
these two effects are visible: at high pj, the global fit has larger resolution tails and worse 
core resolution than the tracker-only fit, while the more selective Tune P performs better than 
tracker-only and sigma switch in both respects. These results are in agreement with previous 
studies and expectations [21. 

6.2.2 Momentum scale from the cosmlcs endpoint method 

The flux of cosmic-ray muons falls steeply as momentum increases. If there is a q/pj bias 
present in the reconstruction at high momentum, the shape of the q/pj spectrum for cosmic- 
ray muons will be distorted, with the location of the minimum shifted from zero; the q/pj bias 
of high-momentum tracks can be extracted from the location of the minimum. The "cosmic 
endpoint" method used to estimate the bias is described in Ref. [|39|| and briefly summarized 
below. 

We perform a binned comparison between the q/pj distributions of the data and of simulated 
events, introducing various artificial biases into the spectrum of simulated events and calculat- 
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ing the as the sum of the squared differences between each bin of the data histogram and the 
corresponding bin of the new simulated histogram, divided by the value in the simulated bin. 
To ensure that the extraction of the bias is not affected by possible data-simulation differences 
in the flux ratio of positive to negative cosmic-ray muons, the c\l 'Pi distribution in the simula- 
tion is normalized to that in the data separately in the q/pj < and q/pj > regions. We then 
fit the distribution of the values as a function of the introduced bias with a polynomial and 
take the location of its minimum as the estimate of the bias in data. The one standard deviation 
statistical uncertainty is computed as the A{q/pj) corresponding to an increase in the function 
value by Ax^ = 1 from its minimum. 

We report results for the tracker-only reconstruction to compare directly with the measure- 
ments performed at intermediate px (Section |6.1| . Using cosmic-ray muons reconstructed with 
Pt > 200GeV/c and passing the selection criteria as described in the previous section, and 
introducing in the simulation a bias k of the form 



q/pj^q/pj + K (9) 

motivated by the arguments used for the choice of the functional forms presented in Eqs. (|2| 
and the bias in data is found to be -0.20 ± 0.12 (stat.) c/TeV. 

To test the robustness of this procedure we generate ensembles of pseudo-experiments using 
simulated cosmic-ray muons for five equally spaced values of k from —0.2 to 0.2 c/TeV. For 
each pseudo-experiment, we repeat the procedure outlined above and calculate a pull as the 
difference between the estimated bias and the input bias divided by the uncertainty. The distri- 
butions of these pulls are consistent with a Gaussian of zero mean and unit width. Considering 
all the pull distributions tested, we assign the largest deviation from zero mean, 0.02 c/TeV, as 
a systematic uncertainty on the estimate of the bias. We also consider the following sources of 
systematic uncertainty and find their impact to be negligible: the resolution model in the sim- 
ulation, the functional form and the fit range used in finding the x^ minimum, and uncertainty 
in the charge ratio of cosmic-ray muons. 



To compare this measurement to that obtained using muons from Z decays (Section 6.1 1, we 
divide the cosmic-ray muon sample into three bins in (p and repeat the above procedure. The 
results using these (p bins and for the whole studied (p range are shown in Table |5j The results 
for the momentum bias at high pj obtained with the endpoint method are compatible with 
those observed at intermediate pj. The accuracy of the measurement is currently statistically 
limited by the available sample of high-px cosmic-ray muons and will improve as more cosmic 
muons are collected. 



7 Background from Cosmic-Ray and Beam-Halo Muons 

7.1 Cosmic-ray muons 

As shown in the previous section and in Ref . IH, cosmic-ray muons are very valuable for study- 
ing the performance of the muon system and muon reconstruction tools. At the same time, they 
are a source of several types of background in physics analyses: 

• A cosmic-ray muon passing close to the interaction point can be reconstructed as a 
collision muon, or as a pair of oppositely charged muons in the upper and lower 
halves of CMS. 

• A muon that is not reconstructed in the tracker (either because it is out-of-time or 
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Table 5: The average c\l pj bias obtained with the cosmic endpoint method in three bins of the 
azimuthal angle cp and for the whole studied (p range. The first uncertainty quoted is statis- 
tical; the second is systematic. Also shown are results obtained with the MuScleFit method 
(Section |6.1 1 by rescaling pj of each of the muons from Z decays according to Eq. ^ and com- 
puting the average difference between the rescaled and the original q/pj values in bins of (p. 
The distribution of muons from Z decays was reweighted to have the same (p distribution as 
that of cosmic-ray muons. The uncertainty quoted is statistical only. 



(p range, rad 


Endpoint method, c/TeV 


MuScleFit method, c/TeV 


-n <^ < -2.1 
-2.1 <(p< -1.05 
-1.05 < ^ < 


-0.14 ±0.09 ±0.02 
-0.31 ± 0.09 ± 0.03 
-0.09 ±0.07 ±0.01 


-0.10 ±0.01 
-0.11 ±0.01 
-0.06 ±0.01 


-TT < (/) < 


-0.20 ±0.12 ±0.02 


-0.10 ±0.01 



passes too far from the interaction point) can still be reconstructed as a standalone 
muon in the muon system and accidentally matched to a tracker track, forming a 
mismeasured global muon. 

• A cosmic-ray muon can deposit energy in the calorimeters but avoid detection in 
the tracking detectors, which would result in, e.g., mismeasured missing transverse 
energy. 

To discriminate between cosmic-ray and collision muons, several quantities can be used, re- 
lated to either properties of the whole event (vertex quality and track multiplicity) or properties 
of a muon candidate (transverse impact parameter \d^y\, timing information, and coUinearity 
of two tracks reconstructed in the upper and lower halves of the detector and associated with 
the same muon). Their performance was studied using two different data samples: a sample of 
muons from collisions with a small (of the order of 1%) contamination from cosmic-ray muons, 
selected by requiring events with a high-quality reconstructed primary vertex, and a sample 
of cosmic-ray muons with a contamination from collision events smaller than 0.1%, selected 
by requiring events without a well-reconstructed primary vertex and with at most two tracks 
reconstructed in the inner tracker. For comparisons with the data, Monte Carlo samples of 
Z — ?► and Drell-Yan dimuon events were used as representative of muons from colli- 

sions for the variables under study. Only events with at least one Global Muon reconstructed 
with Pi > lOGeV/c and passing loose quality criteria were selected. All distributions were 
normalized to the number of collision muons in data. 

A few examples of variables that can be used to distinguish cosmic muons from collision muons 



are shown in Fig. 25 for cosmic-ray muons, muons from collision data, and muons in simulated 



Z and Drell-Yan events. Figure 25 'a) shows the distributions of the transverse impact parame- 



ter calculated relative to the beam-spot position dxy. The d^y distribution for cosmic-ray muons 
is flat, while for collision muons it is strongly peaked at zero. The typical requirement used in 
muon analyses to suppress the cosmic-muon background is |dxy | < 0.2 cm. For illustration, the 
Z+Drell-Yan MC prediction is shown as well, exhibiting an even narrower transverse impact 
parameter peak than the inclusive collision muon sample (cf. Fig. [7]l and demonstrating that a 
tighter and more effective dxy requirement can be applied in specific cases. 

If a cosmic muon passes close enough to the centre of the tracker, it is expected to be recon- 
structed as two tracks ("legs") in the upper and lower halves of CMS with almost equal pj and 
opposite directions. The distinct topology of cosmic muons is well represented by the opening 
angle a, defined as the largest angle between the inner-tracker track of a muon and any other 
track in the event. The distributions of log^Q(7r — a) for cosmic-ray and collision muons are 
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Figure 25: Distributions of variables used for identification of cosmic-ray muons, shown for 
collision and cosmic-ray muon data samples described in the text, and for Z+Drell-Yan MC 
samples: (a) muon transverse impact parameter with respect to the nominal beam-spot po- 
sition; (b) log^Q(7r — a) (see text); (c) muon timing; (d) timing difference between upper and 
lower muon legs. 



shown in Fig. 25 'b). As expected, for most cosmic muons a track opposite in direction could 



be found, which is not the case for collision muons or muons from simulated Z+Drell-Yan 
samples. The excess in the number of collision muons observed at the negative tail of the dis- 
tribution over that predicted is due to a small admixture of cosmic-ray muons present in the 
data. 

Muon timing information can also be used to discriminate between cosmic-ray muons and mu- 
ons from collisions. The timing distributions observed in data and MC samples are shown in 



Fig. 25 'c). The muon timing is defined as the time at which a muon would pass the interac- 
tion point relative to the time of the bunch crossing. It is calculated from timing information 
available in the muon system under the assumption that a particle is moving at the speed of 
light from the centre of the CMS detector outward. Therefore, for collision muons one should 
expect a distribution peaked around zero with a width determined by the resolution of the 
time measurement. For cosmic-ray muons the actual arrival time distribution is flat but the 
reconstructed distribution is also centred around zero. This is due to the trigger and inner- 
tracker reconstruction being most efficient for in-time particles. The additional peak visible at 
negative time is due to the upper reconstructed legs in events where the lower leg is the source 
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of the trigger. In events where both legs are reconstructed, the timing difference between the 
upper and lower legs is a good discriminator between dimuons from collisions and cosmic-ray 



muons, as can be seen in Fig. 25 'd). 



Cosmic-ray muon backgrounds must be removed from all physics analyses using muons. For 
the vast majority of analyses, a simple impact-parameter requirement |dxy| < 0.2 cm, which 
suppresses the contamination from cosmic muons in a sample of Tight Muons with pj > 
10 GeV/c to below 10^**, is sufficient. Some analyses, in particular searches for new phenomena 
involving very energetic muons, need to employ additional cosmic-ray background rejection 
criteria because the proportion of cosmic muons to collision muons increases as the muon mo- 
mentum increases. For example, in the search for extra charged gauge bosons (W) decaying 
to a muon and a neutrino [40], the standard impact-parameter cutoff at 0.2 cm was replaced 
with a tight cutoff of 0.02 cm. In the search for extra neutral gauge bosons (Z') decaying to 
muon pairs ||4T|, the impact-parameter requirement was complemented with a selection on the 
opening angle, (tt — a) > 0.02 rad. As the amount of accumulated data increases, the reach of 
physics analyses is extending to higher values of dimuon mass and muon momentum; since 
the spectrum of muons predicted by the Standard Model falls more rapidly with increasing pj 
than the cosmic-muon spectrum, more sophisticated methods may become necessary to sup- 
press cosmic-ray muon background to the required level in those kinematic regions. A cosmic- 
muon identification algorithm has been developed that combines several quantities sensitive 
to cosmic muons in order to maximize the cosmic-ray muon selection efficiency and minimize 
misidentification of muons from pp collisions as cosmic-ray muons. This algorithm quanti- 
fies the compatibility of any given muon with the hypothesis that it is a cosmic-ray muon; the 
quantities used are the ones described earlier in this section plus a few others characterizing the 
event activity and whether there are hits in the opposite (upper or lower) half of CMS which 
could be associated to a muon track by the dedicated cosmic-muon reconstruction algorithm. 

Table |6] shows the performance of cosmic-ray muon identification with some typical individ- 
ual selections and that of the cosmic-muon identification algorithm. The efficiency is defined 
as the fraction of events identified as cosmic-ray muons by each algorithm in the sample of 
cosmic-ray muons selected from collision data, as discussed above. The misidentification frac- 
tion is the fraction of simulated Z+Drell-Yan events identified as cosmic-ray muons. The first 
row shows the performance of the back-to-back identification, which requires that the opening 
angle a. > {n — 0.1) rad and that \pji — Pti\I ^pii • Pt2 < 0.1. The second row shows the 
effect of applying an impact-parameter requirement: muons with \d^y\ > 0.2 cm are flagged as 
cosmic-ray muons. This requirement is quite efficient in identifying cosmic-ray muons, while 
maintaining a low misidentification probability for collision muons. In the third and fourth 
rows, we show the performance of the cosmic-muon identification algorithm. This algorithm 
allows the analyzer to tune the identification efficiency and misidentification fraction according 
to the needs of the analysis. The tight selection minimizes the misidentification at the cost of 
some loss of efficiency in cosmic-ray muon rejection. The loose selection maximizes cosmic-ray 
muon rejection efficiency at the cost of some higher misidentification. 

7.2 Beam-halo muons 

Accelerator-induced backgrounds ("beam halo") contribute to a variety of physics analyses. 
The CMS detector can identify beam-halo muons that overlap with collision events and might 
otherwise be considered as part of the collision itself. The CSC system has extensive geometri- 
cal coverage in the plane transverse to the beamline and hence can detect halo muons traveling 
parallel to the beamline as well as muons originating from collisions. An algorithm has been 
developed to identify beam-halo muons using information from the CSCs available at both the 
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Table 6: The efficiency of cosmic-ray muon identification and misidentification fractions for 
some relevant individual quantities and the loose and tight versions of the cosmic-muon iden- 
tification algorithm. Uncertainties are statistical only. 





Efficiency (%) 


Misidentification (%) 


Back-to-back 


89.97 ±0.13 


0.153 ± 0.002 


|c?xy| > 0.2 cm 


99.05 ± 0.04 


0.0045 ± 0.0003 


Cosmic ID tight 


97.58 ± 0.07 


0.0001 ± 0.0001 


Cosmic ID loose 


99.52 ± 0.03 


0.153 ± 0.002 



trigger and the reconstruction level. The following three categories of information have proved 
effective: 1) a dedicated Level-1 beam-halo trigger (based on patterns of hits in the CSC system 
consistent with a charged particle traveling parallel to the beamline); 2) CSC trigger primitives 
(patterns of hits within an individual CSC consistent with a traversing charged particle) that 
are early in time relative to the actual pp collision; and 3) a reconstructed standalone-muon 
track with a trajectory parallel to the beamline. Requiring any one of these conditions leads to 
excellent beam-halo muon identification efficiency. We refer to this as Loose Halo ID. Requiring 
any two of these conditions ensures fewer collision muons are spuriously accepted as halo, but 
at the expense of slightly lower efficiency for true halo muons. This is referred to as Tight Halo 
ID. 

The performance of the beam-halo identification algorithms was studied with 2010 collision 
data and Monte Carlo samples. Simulated beam-halo muons that passed through the barrel or 
endcap calorimeters were found to satisfy the requirements of the Loose Halo ID 96% of the 
time, while satisfying the requirements of the Tight Halo ID 65% of the time. To measure these 
efficiencies in collision data we used £^'** -triggered events. Beam-halo muons can induce a 
large missing energy signal by traversing the barrel or endcap calorimeters at constant (p, and 
thus £^^^ -triggered data provide an enriched sample of beam-halo muons. A pure sample of 
beam-halo muons was selected by requiring no reconstructed collision muon in an event and 
£miss ^ 50 opposite in (p to at least one reconstructed hit in the CSCs. The fraction of muons 
in the selected sample identified as beam halo was found to be 89% for the Loose and 73% for 
the Tight Halo ID. Given the number of simplifying assumptions involved in the simulation 
of beam-halo events, the agreement between the data and the simulation is considered to be 
satisfactory. 

The fraction of simulated minimum-bias events in which a particle produced in a proton- 
proton collision was misidentified as a beam-halo muon was ~5 x 10"^ and ?»10"^ for the 
Loose and Tight Halo ID algorithms, respectively. When evaluated from minimum-bias colli- 
sion data, these misidentification probabilities were found to be somewhat larger but still small, 
with values of ~2 x 10^"^ and ~8 x 10^^, respectively. 

The probability for a beam-halo muon crossing the CSCs to be misidentified as a collision muon 
was also evaluated in data and found to be ?»10^^ for Soft Muons with pj < 5GeV/c and 
< 5 X 10^^ for Soft Muons with higher pj and for Tight Muons. Another important number 
is the fraction of collision events contaminated by beam-halo muons. In 2010, it ranged from 
~10^^ to ^10^^, depending on the LHC fill, for muon-triggered events f42]; the fraction could 
be higher by an order of magnitude for hadronic-triggered events, due to the beam-halo muons' 
role in triggering the event via energy deposited in the calorimeters. Several searches for new 
physics in CMS [43l |44| employed the beam-halo identification algorithms described above to 
veto events with beam-halo contamination. 
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The requirement that a muon is an isolated particle in the event, meaning that the energy flow 
in its vicinity is below a certain threshold, can effectively discriminate muons from the decays 
of W and Z bosons from those produced in heavy-flavor decays and hadron decays in flight. 
Three different isolation algorithms have been studied: 

• Tracker relative isolation (I^^l)- This algorithm calculates the scalar sum of the pj of all 
tracker tracks reconstructed in a cone of radius AR = ^^(A^)^ + {Arj)^ < 0.3 centred 
on the muon track direction. The pj of the muon track itself is not included in the 
sum. For the muon to be considered isolated, the ratio of the px sum to the muon 
track px is required to be below a certain threshold. Track directions and values of px 
are computed at the point of closest approach to the nominal centre of the detector. 

• Tracker-plus-calorimeters (combined) relative isolation (Ilg,„i,)- The discriminating vari- 
able is similar to I^^^, but the numerator of the ratio also includes the sum of energies 
measured in ECAL and HCAL towers found within a cone of radius AR < 0.3 cen- 
tred on the muon track direction. The energy deposits associated with the muon 
track itself are not included in the sum. 

• Particle-flow relative isolation (Ipp). The discriminating variable is the sum of the px of 
all charged hadrons, the transverse energies £x of all photons, and £x of all neutral 
hadrons reconstructed by the particle-flow algorithm [21] within a cone of radius 
AR < 0.4 centred on the muon track direction, divided by the muon track pj. 

The optimization of the cone sizes was performed independently for the different algorithms, 
resulting in the different AR values mentioned above. Each of these algorithms has features 
that suit the requirements of different analyses. For example, the tracker-plus-calorimeters 
relative isolation with a threshold of 0.15 is the algorithm chosen for the measurement of the 
W and Z cross sections ^9\. The search for heavy resonances decaying into muon pairs | |4T] [ 
instead used the tracker relative isolation given that high-energy muons are expected to deposit 
a significant amount of energy in the calorimeters. Photons emitted in final-state radiation can 
lead to energy deposited in the ECAL only. An isolation algorithm that uses energy deposits 
in the inner tracker and the HCAL only is not considered in this paper. Analyses involving t 
leptons [45, 46 1 use particle-flow isolation for the muons from the T decays, since particle flow 
is used for the global event reconstruction in these analyses. 

The efficiency of the various isolation algorithms was measured in data, using a sample of 
muons from Z decays. Two methods were used: tag-and-probe and the Lepton Kinematic 



Template (LKT) method. The tag-and-probe method is described in Section 5.1 The LKT 
method [31 is an extension of the random-cone [47J method. It relies on the assumption that 
the kinematics of muons from decays of W or Z bosons produced in the hard parton scatter- 
ing is unrelated to accompanying interactions of the other partons in the colliding protons (the 
underlying event), which are responsible for the energy flow around the muons. An isolation 
variable can then be computed relative to any specific direction in an event with underlying 
event activity similar to that of a signal event. Events containing a Z decaying into a pair of 
muons, with the reconstructed muon tracks and activity associated with them discarded, were 
used as approximations to underlying events. Directions were drawn from template kinematic 
distributions of muons obtained from simulation. As every event can be re-used multiple times, 
the LKT method provides an important cross-check for the results obtained with the tag-and- 
probe technique, which are statistically limited in some kinematic regions. 
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An almost pure sample of Z — )■ events was obtained by selecting events with a pair of 

oppositely charged muons that form an invariant mass in a mass region around the nominal 
Z mass (between 70 and llOGeV/c^ for the tag-and-probe method and within lOGeV/c^ of the 
nominal Z mass for the LKT method) and with each muon satisfying the Tight Muon identifi- 
cation criteria. For the tag-and-probe method, the tag muon was also required to be isolated 
with Ilf^^ < 0.15. A sample of simulated Z — )■ and Drell-Yan dimuon events passing the 

same selection criteria was used for comparisons with the data, with simulated minimum-bias 
collision events overlaid to reproduce the observed pile-up (PU) distribution. 



Figure 26 shows the efficiency of the various isolation algorithms evaluated on muons with 
20 < px < 50 GeV/c from Z decays as a function of the threshold on the corresponding isolation 
variable. Results obtained with the tag-and-probe (for all three isolation algorithms) and the 
LKT (for Il^l and F^^^,) methods are shown for both data and simulation. Efficiencies in data 
and MC simulation are generally in good agreement, with the difference between the two not 
exceeding 1.5%. A single exception is the results from the LKT method for low values of the 
^comb threshold, where efficiencies from MC simulation are lower than those from data by up 
to 4%. The origin of this discrepancy is under study. The ratios of efficiencies obtained on data 



and on simulation are also reported in Fig. 26 and can be considered as scale factors to correct 



MC results. The agreement between the tag-and-probe and LKT results on data is within 1%. 
The efficiency of the various isolation algorithms for muons from Z decays is shown as a func- 



tion of the muon pj and for two threshold values in Fig. 27 for I^^^, Fig. 28 for F^^^, and 



in Fig. 29 for 7pp. In each of the three figures the ratio between data and MC results is also 
reported. The two threshold values studied for each isolation algorithm are those most com- 
monly used in physics analyses and correspond approximately to the end of the rapid rise of 
the efficiency curve and to the beginning of the efficiency plateau shown in Fig. |26] The results 
obtained with the tag-and-probe and LKT methods agree within the statistical uncertainties 
down to the lowest tested muon pj (5GeV/c). The LKT results, which have very small statis- 
tical uncertainties, indicate that for muon pj as low as 5 GeV/ c, the agreement between data 
and MC efficiencies for I^^^ (^comb) within 5% (10%), while for pj greater than 15 GeV/c the 
agreement is within 1%. 

Another important characteristic of the performance of an isolation algorithm is its power to 
reject muons from typical background events. To estimate the rejection power, a sample en- 
hanced with events from QCD processes is used. These events are selected by requiring the 
following: one offline reconstructed muon with 20 < px < 50 GeV/ c satisfying the Tight Muon 
requirements; no other Tight Muon reconstructed in the event; at least one particle-flow jet 
reconstructed with £x greater than 30 GeV; the azimuthal angle A(p^^ j-miss between the muon 
and the missing transverse energy directions smaller than 1.5 rad; and the transverse mass. 



defined as Mj = ^2pjE™''^^{l — cos(A(/)^^ e^'^^O)' smaller than 20GeV/c^. The last two require- 
ments reduce the contamination from W events significantly. A sample of simulated events 
from QCD processes, W decays, Z — > , Drell-Yan, and ti production, mixed according to 
the corresponding cross sections, is used to study the expected composition of the data sample 
resulting from the above selection and for comparing results obtained from data and MC simu- 
lations that include pile-up observed in data. Before isolation requirements are applied, events 
from QCD processes represent 99.6% of the MC sample. The next most important contributions 
come from W (0.3%) and Z+Drell-Yan events (0.1%). 



Figure 30 shows the efficiency of the various isolation algorithms for muons from Z decays 
versus the efficiency of the same algorithm for muons from the QCD-enhanced dataset. In 
both cases muons are required to have 20 < px < 50 GeV/c. Results obtained with the tag-and- 
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Figure 26: Left: efficiencies of various isolation algorithms for muons with 20 < px < 50GeV/c 
from Z decays as a function of the isolation threshold. Results are shown for both data and 
simulation using the tag-and-probe ("T&P") and Lepton Kinematic Template ("LKT") meth- 
ods; the LKT method is not used for the particle-flow algorithm. Right: data to simulation 
ratios. Plots are shown for tracker relative (I^^l, top), tracker-plus-calorimeters relative (Ilf^^^, 
middle), and particle-flow relative (7pp, bottom) isolation algorithms. The MC samples include 
simulation of pile-up events. 
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Figure 27: Left: efficiency of tracker relative isolation I^^^ for muons from Z decays as a func- 
tion of muon pj, measured with data using the tag-and-probe ("T&P") method (in bins of [5, 
10], [10, 20], [20, 50], and [50, 100] GeV/c) and the LKT method (finer binning). Results corre- 
sponding to the threshold values of 0.05 and 0.1 are shown. Right: ratio between data and MC 
efficiencies. The MC samples include simulation of pile-up events. 
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Figure 28: Left: efficiency of tracker-plus-calorimeters relative isolation Ilf^^ for muons from Z 
decays as a function of muon pj, measured with data using the tag-and-probe ("T&P") method 
(in bins of [5, 10], [10, 20], [20, 50], and [50, 100] GeV/c) and the LKT method (finer binning). 
Results corresponding to the threshold values of 0.10 and 0.15 are shown. Right: ratio between 
data and MC efficiencies. The MC samples include simulation of pile-up events. 

probe method are reported; the LKT results agree with the tag-and-probe results within 1% for 
both data and simulation, and are shown only for the J^o^j^b algorithm for comparison. As can 
be seen in the figure, the best performance for the considered datasets and px range is given by 
7pp, followed closely by Ilf^^^- For a given signal efficiency, the background rejection evaluated 
using the MC sample (not shown) is slightly better than that in the data, by up to 0.5% for I^^^ 
and up to 1% for Ilf^i^ and 7pp. It should be noted that the reported background efficiency 
includes the absolute 0.4% contamination from truly isolated muons from W/Z decays in the 
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Figure 29: Left: efficiency of particle-flow relative isolation 7pp for muons from Z decays as a 
function of muon pj, measured with data using the tag-and-probe ("T&P") method. Results 
corresponding to the threshold values of 0.12 and 0.20 are shown. Right: ratio between data 
and MC efficiencies. The MC samples include simulation of pile-up events. 



background sample, as described above. 
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Background Efficiency 

Figure 30: Isolation efficiency for muons from Z decays versus isolation efficiency for muons in 
the QCD-enhanced dataset described in the text, for tracker relative, tracker-plus-calorimeters 
relative, and particle-flow relative isolation algorithms. Muons are required to have pj in the 
range between 20 and 50 GeV/c. The background rejection is limited by the 0.4% contamination 
from truly isolated muons. 

The efficiency of any isolation algorithm features a dependence on the muon pseudorapidity, 
and this dependence is expected to become more pronounced as the number of pile-up col- 
lisions increases. Figure 31 shows the efficiency for muons with 20 < px < 50GeV/c from 



Z decays observed in data compared with the results obtained using MC samples having the 
same distribution of the number of primary vertices as in data, with an average of 2.4. The 
results are obtained with the LKT method using the Ilf^^ algorithm with a threshold of 0.15. 
The particular t] dependence reflects the distribution of transverse energy flow expected for 
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Figure 31: Muon isolation efficiency versus rj. Efficiency for muons with pj between 20 and 
50 GeV/c from Z decays observed in 2010 data, in which the average number of reconstructed 
primary vertices (PV) is 2.4, is compared with the results obtained using MC samples with the 
same average number of primary vertices. The results were obtained with the LKT method 
using the I^f^^^ algorithm with a threshold of 0.15. 

minimum-bias events. Methods to mitigate the impact of pile-up on the performance of iso- 
lation algorithms have been developed for higher luminosity running. One such technique 
is based on the measurement, event-by-event, of the average transverse momentum per unit 
area p added to the event by minimum-bias pile-up collisions ||48ll . Another technique uses 
reconstructed tracks and primary vertices to compute a correction factor j6 to be applied to all 
or part of the numerator of the isolation variables described above. In both cases, the p and /3 
variables allow the energy in the isolation cone due to particles produced in pile-up collisions 
to be estimated. 

An additional advantage of the LKT method is that it can be used to estimate the isolation 
efficiency for muons produced in any signal process. The basic assumption is made that the 
kinematics of a muon produced in the leptonic decay of a W, Z, or a heavier new particle is 
unrelated to the activity of the rest of the event. Under this assumption, the isolation efficiency 
for muons from such a signal can be inferred from the corresponding efficiency measured using 
muons from Z, by extrapolating from underlying event activity of Z events to that of the signal. 
The number of ECAL and HCAL calorimeter towers with transverse energy above threshold 
can be used as a measure of the event activity. The power of this technique can be appreciated 
in Fig. 32 where the isolation efficiency for prompt muons from simulated tt events forced 
to decay leptonically is plotted versus the number of calorimeter towers with < 2.7 and 
Ej > 1 GeV. The isolation efficiency for muons from Z decays (data and MC simulation) is 
shown superimposed. The efficiencies shown are those of the I^f^}^ algorithm computed by 
the LKT method. Compared to tt, Z events are characterized by relatively low occupancies. In 
the region where the Z and tt efficiencies overlap, however, the differences between them do 
not exceed 2%, and both sets of points tend to align along an approximately straight line. The 
parameters of this line can be obtained by fitting the points from the Z dataset; the parameters 
are then used to predict, with an accuracy of a few percent, the isolation efficiency for muons 
in events with large underlying activity, which cannot be directly measured with the tag-and- 
probe technique. The pile-up contribution can then be accounted for and corrected by means 
of the techniques using the p and /3 variables described above. 
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Figure 32: Muon isolation efficiency versus number of calorimeter towers in the region of | // 1 < 
2.7 having transverse energy larger than 1 GeV. Results for muons with 20 < px < 50GeV/c 
from Z decays obtained in data and in simulation are compared with the efficiency for prompt 
muons in a simulated sample of top pairs decaying leptonically. All results were obtained with 
the LKT method using the I^^^^, algorithm with a threshold of 0.15. The MC samples include 
simulation of pile-up events. 

9 Muon Trigger 

The CMS trigger system consists of two basic stages: the hardware-based Level-1 trigger 
and the software-based high-level trigger (HLT) [6J. 

The Level-1 muon trigger uses signals from all three CMS muon detector systems: DT, CSC, and 
RFC. It has a latency of 3.2 jis and reduces the rate of inclusive muon candidate events read-out 
from detector front-end electronics to a few kHz by applying selections on the estimated muon 
pj and quality. 

In the muon HLT, first a Level-1 trigger object is used as a seed to reconstruct a standalone- 
muon track in the muon system, leading to an improved px estimate. At this point, px thresh- 
old filters are applied to the standalone (also called Level-2) muon. Then seeds in the inner 
tracker are generated in the region around the extrapolated Level-2 muon, and tracker tracks 
are reconstructed. If a successful match is made between a tracker track and the Level-2 muon, 
a global fit combining tracker and muon hits is performed, yielding a Level-3 muon track on 
which the final pj requirements are applied. In this way, the rate of recorded inclusive muon 
events is reduced to a few tens of Hz. The average processing time of the HLT reconstruction 
is about 50 ms. 

In this section we report on trigger efficiency measurements and rejection rates, and compare 
them to predictions from the simulation. We study the trigger efficiency using two comple- 
mentary methods: one using tag-and-probe and one using single muons reconstructed offline. 
The efficiencies are shown separately for the barrel and for the overlap-endcap regions since 
differences in the trigger response have been observed between them. 

9.1 Trigger efficiency using the tag-and-probe method on dimuon resonances 



In this subsection, we report the measurements of the trigger efficiency for prompt muons 
performed by applying the tag-and-probe method described in Section 5.1 to muons from the 
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decays of] /ip and Z resonances. With this method, the trigger efficiency can be measured with 
respect to any offline muon selection. In the following, we use Soft Muons and Tight Muons 
as probes to check how efficient the trigger is in selecting muons that can be reconstructed 
offline. In the case of the Z, the probe muons are also required to be isolated by requiring I^^^^ 
to be smaller than 0.15. Measured efficiencies are compared to those evaluated using the MC 



samples described in Section 5.1 



The muon trigger efficiency was studied separately in two regions of muon transverse momen- 
tum: below and above 20GeV/c. In the region of pj > 20GeV/c, the efficiency was measured 
using a sample of Z — > ji^ji^ events collected with single-muon triggers. In the region of 
Pt < 20GeV/c, the ]/ip V-^l^^ events collected by the specialized muon-plus-track triggers 
described in Section |2] were used. Several instances of the muon-plus-track trigger have been 
deployed, with different thresholds on the px of the tracker track; those with the lower thresh- 
olds became prescaled as the instantaneous luminosity increased during 2010. As a result, the 
amount of data used for the measurements presented here varies with the px region: it corre- 
sponds to an integrated luminosity of 0.9 pb^^ for px < 3 GeV/c, 0.6 pb^^ for 3 < px < 5 GeV/c, 
and 5.3 pb^^ for 5 < px < 20GeV/c. The sample of Z — > ji^^^ events used for the measure- 
ments in the range of px above 20 GeV/c corresponds to an integrated luminosity of 31 pb^^. 

To evaluate the trigger efficiency, trigger objects must be matched to the muons reconstructed 
offline. Level-1 muon-trigger candidates are matched to offline muons by position, extrapo- 
lating a muon tracker track to the muon system. The HLT muons are matched to the muons 
reconstructed offline by direction at the vertex. The HLT-only efficiencies are computed match- 
ing the probe with the Level-1 candidate and requiring it to be also matched with the HLT 
candidate. The combined efficiencies of Level-1 and HLT are obtained simply requiring the 
probe to be matched with the HLT candidate. 

9.1 .1 Trigger efficiency for Soft Muons 

Trigger efficiencies for Soft Muons in the range of px < 20GeV/c were measured using muons 
from the decays of ]/\p particles. 



The single-muon trigger efficiencies for Soft Muons are shown in Fig. 33 as a function of muon 
Px, separately for the Level-1 trigger with px threshold at 3 GeV/ c, the HLT with px threshold 
at 5 GeV/ c, and for the full Level-l-HLT online selection. The efficiencies at the plateau, for 
muons with px in the range of 9 < px < 20 GeV/ c, are reported as a function of muon rj in 
Fig. 34 and summarized in Table [7| 



All efficiency curves show a sharp "turn-on" at the trigger px threshold, well described by the 
simulation. As expected, the turn-on is sharper for the HLT, due to an improved px resolution. 
The plateau efficiencies are about 99% for the Level-1 trigger in the barrel region and for the 
HLT in the whole studied pseudorapidity range. The Level-1 efficiency in the overlap-endcap 
region is slightly lower, at about 95%. Most of the 5% efficiency loss is due to stringent quality 
criteria used in the selection logic of Level-1 muon triggers during 2010 data taking; these cri- 
teria were further optimized during the 2010-11 winter technical stop of the LHC. Measured 
efficiencies are generally in good agreement with those expected, with the ratios between the 
two (see Table [tJi being within 1% of unity. The only exception is a slightly larger, about 2%, dif- 
ference between data and simulation in the overlap-endcap region due to a few non-operational 
CSCs not accounted for in the simulation (see Section|5.1|. 



The very forward region, 2.1 < |//| < 2.4, is characterized by higher rates of low-px muons 
and poorer momentum resolution. At the nominal LHC luminosity, this region is intended to 
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Figure 33: Single-muon trigger efficiencies for Soft Muons as a function of the Soft Muon pj 
in the barrel (top) and the overlap-endcap (bottom) regions: the efficiency of the Level-1 trig- 
ger with pt threshold at 3GeV/c (left), the efficiency of the HLT with pj threshold of 5GeV/c 
with respect to Level-1 (middle), and the combined efficiency of Level-1 and HLT (right). The 
efficiencies obtained using J/i/? — ?► events (points with error bars) are compared with 

predictions from the MC simulation (ztlc bands); the uncertainties are statistical only. 



be used to improve trigger efficiency for events with multiple muons but not to be included 
in triggering on single muons. In 2010, however, the luminosity and hence the rates of low- 
Pj muons were sufficiently low to allow the single-muon trigger to be extended to include 
the entire acceptance of muon detectors, up to \rj\ = 2.4. In this region of \rj\ = 2.1-2 A, the 3:1 
ganging of 48 CSC cathode strips into 16 readout channels in the first muon station (49] leads to 
an ambiguity in pj assignment at Level-1 trigger. The algorithm to resolve this ambiguity was 
configured with the goal to reduce the rate of muons with overestimated pj to a minimum. 
As a consequence, the efficiency of the Level-1 triggers in the very forward region was high 
for triggers with low pj thresholds (about 95% for the Level-1 trigger with pj threshold at 
3GeV/c, see Fig. |34]|, but substantially lower for triggers with thresholds above 3GeV/c (see 
next section). 

9.1 .2 Trigger efficiency for Tight Muons 

Trigger efficiencies for Tight Muons were measured by applying the tag-and-probe method to 
]/xp — ?► ji^ji^ events in the region of pj below 20 GeV/c and Z — )■ events in the region of 

Pt above 20 GeV/c. 

Using this combination of these two resonances, we can study efficiencies of the Level-1 trigger 
with Pt threshold at 7GeV/c and of two high-level triggers, with thresholds at 9 and 15 GeV/c. 
These efficiencies, together with the combined efficiency of the Level-1 and HLT, are shown in 



Fig. 35 as a function of muon pT. The efficiencies for muons with pT > 20GeV/c evaluated 
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Figure 34: Single-muon trigger efficiencies for Soft Muons as a function of the Soft Muon t], for 
muons with pj in the range of 9 < px < 20 GeV/ c: the efficiency of the Level-1 trigger with pj 
threshold at 3 GeV/ c (left), the efficiency of the HLT with pj threshold of 5 GeV/ c with respect 
to Level-1 (middle), and the combined efficiency of Level-1 and HLT (right). The efficiencies 
obtained using ]/xp — > j^^H^ events (points with error bars) are compared with predictions 
from the MC simulation (±1(7 bands); the uncertainties are statistical only. 



Table 7: Level-1, HLT, and overall single-muon trigger efficiencies for Soft Muons in the ef- 
ficiency plateau (9 < px < 20GeV/c) for different pseudorapidity regions. The first column 
shows efficiencies measured from data; the second column shows the ratio between the mea- 
surements in data and simulation. The uncertainties are statistical only. 



Trigger Level 


Region 


Tag-and-Probe J/t/j 
Eff. [%] Data/MC 


Level-1 \r]\ < 2.1 
\t]\ < 0.9 
0.9 < \f]\ < 2.1 


97.1 ± 0.2 0.990 ± 0.002 

99.2 ± 0.1 0.995 ± 0.001 
94.5 ± 0.3 0.978 ± 0.004 


HLT l'7l<2.1 
\ri\ < 0.9 
0.9 < \r]\ < 2.1 


99.1 ± 0.2 0.995 ± 0.002 

99.2 ± 0.2 0.993 ± 0.002 
99.0 ± 0.3 0.997 ± 0.004 


Level-l+HLT \?]\ < 2.1 
\ri\ < 0.9 
0.9 < \r]\ < 2.1 


96.2 ± 0.2 0.985 ± 0.003 

98.5 ± 0.2 0.989 ± 0.002 

93.6 ± 0.5 0.975 ± 0.005 



events are shown as a function of muon rj in Fig. 36 and are summarized in 



using Z 

the first two columns of Tabled 

Like the trigger efficiencies for Soft Muons, the efficiency curves for Tight Muons show a rapid 
turn-on of efficiencies near the applied threshold: for Level-1 threshold at pi = 7GeV/c and 
HLT threshold at px = 9GeV/c, the plateau is reached at px ~ lOGeV/c. The turn-on region 
is well reproduced by the simulation. The combined Level-1 and HLT efficiency at the plateau 
is about 95% in the region of < 0.9 and about 90% in 0.9 < \r]\ < 2.1. Most of the effi- 
ciency loss occurs at Level-1 and, in particular, in the overlap-endcap region for the reasons 
described in Section 9.1.1 In the very forward region, 2.1 < \f]\ < 2.4, the single-muon trigger 
efficiency is about 45%. Most of the losses are due to the underestimation of px by the Level-1 



px-assignment algorithm described in Section 9.1.1 this algorithm has been modified during 
the 2010-11 winter technical stop of the LHC to yield higher efficiency for muons with px larger 
than 20 GeV/c. The difference between the measured and predicted average plateau efficiencies 
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Figure 35: Single-muon trigger efficiencies for Tight Muons as a function of the Tight Muon pj 
in the barrel (top) and the overlap-endcap (bottom) regions. The measurements are done with 
the tag-and-probe method, using ]/tp H^H^ events for pj below 20GeV/c and Z — > 
events above. The efficiencies are shown for the following triggers: Level-1 with pj > 7GeV/c 
threshold (left), HLT with pj threshold at 9 GeV/c for pj below 20 GeV/c and at 15 GeV/c above 
(centre), and the combination of the above Level-1 and HLT triggers (right). The efficiencies in 
data (points with error bars) are compared with predictions from the simulation (±1(7 bands); 
the uncertainties are statistical only. 



(see Table|8]l is within 0.5% for the HLT and of the order of 2-3% for the Level-1. 

To estimate the effect of multiple interactions on the muon trigger performance, the measure- 
ment of the efficiency at the plateau is performed as a function of the number of reconstructed 
primary vertices in the event. Just as for muon identification, no loss in efficiency has been 
observed for events containing up to six reconstructed primary vertices. 

9.1 .3 Systematic uncertainties 

The studies of systematic uncertainties in the data-to-simulation ratios of trigger efficiencies 
closely followed those performed for the muon reconstruction and identification efficiencies 
(Section [5. 1.3| and focused on lineshape modelling and background subtraction. 

The efficiencies obtained by applying the tag-and-probe method to the simulated samples of 
muons are compared with the "true" efficiencies computed by simple counting of the passing 
and failing probes in simulated ^/tp — )■ fi^ji^ and Z — )■ fi^ji^ events. The results are in good 
agreement: the difference in the efficiencies is smaller than 0.1% for muons in both samples, 
well within the statistical uncertainties of the measurements. 

For ]/tp — )■ fi^ji^ events in data, the efficiencies are recomputed using a simple Gaussian in- 
stead of a Crystal Ball function [28J to model the resonance and with a quadratic polynomial 
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Figure 36: Single-muon trigger efficiencies for Tight Muons with pj > 20 GeV/ c as a function 
of the Tight Muon rj: the efficiency of the Level-1 trigger with pj threshold at 7GeV/c (left), 
the efficiency of the HLT with pj threshold of 15GeV/c with respect to Level-1 (middle), and 
the combined efficiency of Level-1 and HLT (right). The efficiencies obtained using Z — )■ 
events (points with error bars) are compared with predictions from the MC simulation (±1(7" 
bands); the uncertainties are statistical only. 



Table 8: Summary of Level-1, HLT, and overall trigger efficiencies for Tight Muons at the effi- 
ciency plateau in different pseudorapidity regions. The results obtained using muons from Z 
decays and single isolated muons collected with jet triggers are averaged over the pj range of 
20-100 GeV/c; the results from single muons contained in b-tagged jets are reported without an 
upper Pj limit. For each method, the first column shows efficiencies measured from data; the 
second column shows the ratios between the estimates in data and in simulation. The quoted 
uncertainties are purely statistical. 



Trigger Level 
Region 


Z ji^ji 
Eff. [%] Data/MC 


Single isolated }i 
Eff. [%] Data/MC 


Single in b jets 
Eff. [%] Data/MC 


Level-1 

\n\ < 2-1 

\t]\ < 0.9 
0.9 < ?/ < 2.1 


94.1±0.2 0.976±0.002 
95.9±0.2 0.981±0.002 
92.3±0.3 0.971±0.004 


92.9±0.8 0.966±0.011 
94.7±1.0 0.971±0.013 
91.1±1.2 0.961±0.017 


94.1 ±0.3 0.975±0.003 
95.4±0.3 0.978±0.004 
92.2±0.5 0.971±0.006 


HLT 

<2.1 
|?/| < 0.9 
0.9 < // < 2.1 


98.2±0.1 0.996±0.001 
99.0±0.1 0.996±0.001 
97.4±0.2 0.995±0.002 


97.8±0.5 0.993±0.007 
98.6±0.5 0.996±0.008 
97.1±0.8 0.991±0.011 


94.4±0.3 0.977±0.003 
96.9±0.3 0.988±0.003 
90.5±0.6 0.959±0.007 


Level-l+HLT 
\f]\<2.1 
< 0.9 
0.9 < \t]\ < 2.1 


92.5±0.3 0.972±0.003 
95.0±0.3 0.978±0.003 
89.9±0.4 0.966±0.004 


90.8±0.9 0.959±0.013 
93.3±1.1 0.967±0.015 
88.5±1.4 0.952±0.020 


88.9±0.4 0.953±0.004 
92.5±0.4 0.967±0.005 
83.5±0.8 0.931±0.008 



instead of an exponential to model the background. Like the results reported in Section |5| the 
differences in the efficiencies resulting from this variation in assumed signal shape are smaller 
than 0.1%; the changes in efficiencies due to a different background parametrization are less 
than 0.2%. The same test has been made for Z — > ji^ji^ events: using a quadratic polynomial 
instead of an exponential to model the background results in the efficiency changes of 0.05%. 
In all cases, the difference between the two results gives an estimate of possible systematic 
uncertainty due to the modelling of signal and backgrounds. 

Finally, for the efficiencies calculated using ]/^ events, the uncertainty due to residual corre- 
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lation effects between the two muons has been studied by changing the separation criteria as 
described in Section l5.1.3[ The effect on the measured ratio of efficiencies in data and in simu- 
lation is about 0.3% in the barrel and less than 0.1% in the endcaps. 



9.1 .4 Efficiency of online muon isolation requirements 

To further reduce the trigger rate, isolation criteria can be applied at the HLT. The isolation 
requirements are based on the calorimeter information for the Level-2 muons and on the pixel 
tracks for the Level-3 muons ||6l. 

The efficiency of online isolation requirements was evaluated by applying the tag-and-probe 
technique to muons from Z decays. To obtain the trigger efficiency for muons that are com- 
monly used in physics analyses, the probe is required to pass the Tight Muon selection. Further- 
more, to determine the additional effect of online isolation criteria, we require that the probe 
muon match the trigger object that passed a non-isolated muon trigger of the same threshold. 



Figure 37 shows the efficiency of the online isolation selection as a function of the threshold 
applied to each of the three offline isolation variables described in Section |8] the tracker ab- 
solute isolation (the numerator of the tracker-plus-calorimeters relative isolation (Ilf-^^^), 
and the particle-flow relative isolation (/pp). The probes are Tight Muons with \f]\ < 2.1 and 
Pj > 20 GeV/c, matched with the muon trigger candidates that passed the Level-1 trigger with 
a Pj threshold of 7GeV/c and HLT with a pj threshold of 9GeV/c. We then ask whether the 
matched trigger object passes the combination of the Level-1 and HLT triggers with the same 
pj thresholds, but also including the isolation criteria. 
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Figure 37: Efficiency of single-muon trigger including isolation requirements for Tight Muons 
matched with the trigger objects passing trigger without isolation. Efficiencies measured in 
data (circles) are compared with the predictions of the MC simulation (triangles) as a function of 
the thresholds applied to the offline isolation variables: tracker absolute isolation (left), tracker- 
plus-calorimeters relative isolation (middle), and the particle-flow relative isolation (right). 



The efficiency as a function of the threshold on the tracker-based isolation shows a drop around 
1 GeV/ c. This is a direct consequence of the thresholds implemented in the pixel-track-based 
isolation at Level-3. For the typical values of the thresholds applied offline (see Section |8]l, 
an additional efficiency loss due to online isolation requirements is of the order of 4-5%. The 
measured efficiency distributions are well described by the simulation, with differences of the 
order of 1% or less over the wide range of threshold values. 



9.2 Muon trigger efficiencies from jet-triggered samples 
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9.2 Muon trigger efficiencies from jet-triggered samples 

While the tag-and-probe technique using muons from j/tporZ allows one to determine trigger 
efficiencies for prompt signal muons with sufficiently small bias, the momentum range over 
which the efficiency can be probed is restricted to the momentum range covered by muons 
from the decays of these resonances. 

As an alternative to the tag-and-probe method, it is useful to study the inclusive trigger ef- 
ficiency for single muons reconstructed offline, in a sample collected using a trigger not in- 
volving the muon and tracker systems. In this section, we report trigger efficiencies obtained 
from 2010 data samples recorded with jet triggers, which use only energy measurements in 
the calorimeters. We use data collected using two single-jet triggers, with thresholds on un- 
corrected jet Pt of 70 and lOOGeV/c, as well as events in which the sum of all uncorrected 
transverse energy Hj from jets is larger than 100 GeV. To reduce background, the probes are 
required to be Tight Muons. They are matched to the HLT muon using the track direction at 
the vertex. Mismatches are reduced by requiring that only one reconstructed muon be present 
in the event. 

Backgrounds from pion and kaon decays, as well as remaining hadron punch-through, lower 
the estimated efficiency. From Table |2] we expect these backgrounds to be at the level of about 
10% for Tight Muons with pj > 20GeV/c. A higher rejection of these backgrounds can be 
achieved by applying two alternative and independent additional selections: 

• Isolation, which selects mainly muons from W decays and suppresses the heavy- 
flavour contributions that are typically non-isolated. 

• B tagging, which selects muons from semileptonic heavy-flavour decays. Jets are 
reconstructed by the particle-flow algorithm Il2ll| , and a b-tagging criterion 1150 1 is 
imposed on the jet associated with the muon by requiring a secondary vertex formed 
from at least two high-quality tracks. 



Trigger efficiencies obtained using isolated muons are shown in Fig. 38 as a function of the 
reconstructed muon pj, both for data and for simulated W — >• decays. To select isolated 
muons, the tracker-plus-calorimeters relative isolation is required to be 11^^^, < 0.15. The pj 
spectrum of muons in the sample thus obtained extends beyond 300 GeV/c, thus probing a 
significantly wider momentum range than the tag-and-probe method. The efficiency beyond 
Pt > 20GeV/c is consistent with being flat. The average plateau efficiencies are reported in 
the third and fourth columns of Table |8] They can be directly compared to the results from the 
tag-and-probe method at the Z — > resonance, which effectively uses isolated muons as 

well. The results from the two methods are in good agreement. 

The trigger efficiencies for Tight Muons in jets, obtained from the b-tagging approach, are 



shown in Fig. 39 for data and for a muon-enriched sample of minimum-bias events generated 
using PYTHIA. The simulation includes the same jet-trigger requirements as used in the data. 
The shape of the efficiency curves is well reproduced by the simulation. In the plateau above 
20GeV/c the data and simulation agree to about the 5% level. The plateau trigger efficiency 



as a function of the muon pseudorapidity is shown in Fig. 40 and summarized in the last two 
columns of Table |8l 

The systematic uncertainty due to possible residual background from hadron punch-through 
and muons from light flavours surviving the b-tagging selection was studied by varying the 
selection criteria. The results do not change significantly if 

• before applying the b-tagging condition, the offline muon requirements are further 
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Figure 38: Single-muon trigger efficiencies for Tight Muons passing the offline isolation se- 
lection as a function of the muon pj, in the barrel (|7/| < 0.9, left) and the overlap-endcap 
(0.9 < \ri\ < 2.1, right) regions. In both regions the efficiencies for the following triggers are 
shown: the Level-1 trigger with pj threshold at 7GeV/c (top), the HLT with pj threshold at 
15 GeV/c with respect to the Level-1 (middle), and the combination of above Level-1 and HLT 
triggers (bottom). The efficiencies in data (points with error bars) are compared with predic- 
tions from the simulation (±1(7 bands); the uncertainties are statistical only. 



tightened by imposing more stringent selection criteria on the of the tracker-only 
and the global track fits, and by increasing the minimum required number of pixel 
and strip tracker hits; 

• the b-tagging condition is tightened, by requiring at least three instead of the default 
two high-quality tracks associated with the secondary vertex; 

• only those events are used in which the muon track is matched to the secondary 
vertex by the b-tagging algorithm. 

Performing the same checks on the sample of simulated events leads to similar results. In ad- 
dition, we verified in simulation that the efficiency determined for muons from heavy-flavour 
decays agrees with that determined by applying the b tagging to the full Monte Carlo sam- 
ple including backgrounds. This demonstrates that b tagging is very effective in suppress- 
ing the background. However, due to pre-selection requirements applied at the generator 
level, the simulated sample of muon-enriched minimum-bias events used in this study lacks 
punch-through and light-hadron decays occurring downstream of the ECAL, which results in 
a smaller impact of the b tagging on the muon selection in simulation than in data. Therefore 
we consider the difference between the results obtained with and without the b tagging, which 
is smaller than 1% both in the barrel and endcap, to be a conservative estimate of the systematic 
uncertainty due to background. 

The trigger efficiencies determined for muons in b jets are lower than those obtained for isolated 
muons in similar ranges of muon transverse momentum. The dependence of trigger efficiency 
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Figure 39: Single-muon trigger efficiencies for Tight Muons contained in b-tagged jets as a 
function of the muon pj, in the barrel (top) and the overlap-endcap (bottom) regions. The effi- 
ciencies are shown for the following triggers: the Level-1 trigger with pj threshold at 7GeV/c 
(left), the HLT with pj threshold at 15 GeV/c with respect to the Level-1 (middle), and the com- 
bination of above Level-1 and HLT triggers (right). The efficiencies in data (points with error 
bars) are compared with predictions from the simulation {±la bands); the uncertainties are 
statistical only. 
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Figure 40: Single-muon trigger efficiencies for Tight Muons contained in b-tagged jets as a 
function of the muon pseudorapidity, for muons with pj > 20GeV/c: the efficiency of the 
Level-1 trigger with pj threshold at 7 GeV/ c (left), the efficiency of HLT with pj threshold at 
15 GeV/c with respect to the Level-1 (middle), and the combined efficiency of Level-1 and HLT 
(right). The efficiencies in data (points with error bars) are compared with predictions from the 
simulation (ztlc bands); the uncertainties are statistical only. 



on isolation is demonstrated in Fig. 



41 



which shows Level-1 and HLT efficiencies for muons 
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in a jet- triggered sample versus the tracker-plus-calorimeter relative isolation J^pJ^j,. The thin 
vertical lines at 0.15 show the typical value defining the isolated muon selection. Without the 
b tagging, a slight dependence of the Level-1 efficiency on isolation can be seen, due to the 
residual background contamination; the efficiency becomes flat once the b tagging is applied. 
The HLT efficiency, on the other hand, has a clear dependence on the muon isolation, even 
after b tagging. The downward trend has a steeper slope without the b-tagging condition 
due to the larger background. Similar conclusions are obtained by studying other isolation 
variables, e.g., the track multiplicity within an angular cone around the muon direction. A 
worse performance of the single-muon trigger at HLT for muons within high-multiplicity jets 
is not unexpected. The precision in the matching of the Level-2 muon candidate from the muon 
chambers with the tracker track from the inner tracker is limited by time constraints dictated by 
the online operation. Hence, with higher track multiplicity in the tracker matching region the 
probability of wrong association increases. Nevertheless, the HLT efficiency remains greater 
than 90% even in the most unfavourable environment of high multiplicities, corresponding to 
the highest values of J^^n^i, 
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Figure 41: Single-muon trigger efficiencies for Tight Muons with pj > 20GeV/c and \rj\ < 2.1, 
with and without the further b-tagging selection, as a function of the muon isolation variable 
^comb- Level-1 trigger (left), and the HLT with respect to the Level-1 (right). 



The ratios of muon trigger efficiencies derived from jet-triggered events in data and in simu- 
lation are shown in the last column of Table |8] Since muons in jet-triggered events tend to be 
less isolated than those in muon-triggered events, these corrections are not directly applicable 
to the muon-triggered events and must themselves be corrected for the difference in isolation 
in the two samples. This is done by reweighting the jet-triggered events to match the distri- 
bution of a given isolation quantity in the muon-triggered events, and then recalculating the 
corrections. After the reweighting is done, the data-to-simulation efficiency ratios are closer to 
unity than those before reweighting by 0.01 in the region of |?/ 1 < 0.9 and by 0.03 in the region 
of 0.9 < < 2.1. These final ratios are in agreement with those obtained with the other two 
methods, which are shown in Table [s] Reweighting the events using either the I^^^^, distribu- 
tion or the distribution of the number of tracks in the vicinity of the muon gives statistically 
compatible results. 

The possible effect of event pile-up on the trigger efficiencies was examined by performing the 
efficiency measurements as a function of the number of good reconstructed primary vertices in 
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the event. No significant dependence, within the statistical uncertainty of about 2%, was found 
in the available range of one to seven reconstructed primary vertices. 

9.3 Trigger rejection rates 

The main task of the trigger is to reduce the rate of events to be recorded, while keeping high 
efficiency for the physics-signal events. The two previous sections focused on the trigger effi- 
ciency; in this section we show what fraction of minimum-bias events is rejected as a function 
of muon trigger pj threshold at Level-1, Level-2, and Level-3, and compare to the prediction 
from simulation of minimum-bias events. 

The rates of accepted muon triggers were evaluated on a data sample collected using the Level- 
1 muon trigger with a pj threshold of 7 GeV/ c and corresponding to an integrated luminosity 



of 26.5 nb ^. The resulting cross sections are shown in Fig. 42 At the peak instantaneous 



luminosity of about 2 x 10 cm^ Hz, achieved at the end of the LHC operation in 2010, the 
rate of Level-1 muon trigger with the px threshold of 7GeV/c was about 2 kHz and the rate of 
the HLT trigger with the pj threshold of 15 GeV/c was slightly below 20 Hz. 
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Figure 42: The accepted cross section of events as a function of muon trigger pj threshold for 
the actual Level-1, Level-2, and Level-3 muon trigger objects processed online in data, com- 
pared to the emulated Level-1 and HLT trigger in simulation. Points corresponding to thresh- 
olds of prescaled Level-3 triggers are not plotted. 



To reproduce the trigger rate, the simulation must reproduce the correct sample composition 
and the correct trigger efficiency, and also must correctly describe the resolution for muons 
from different sources including the resolution tails. The level of agreement between the Level- 
1 and HLT trigger performance in data and the Monte Carlo simulation is generally fairly 
good, with the simulation slightly underestimating the rates at higher thresholds. These re- 
sults demonstrate that the trigger simulation is a useful tool to predict muon trigger rates and 
confirm that the performance of the CMS muon trigger system is close to the design expecta- 
tion. 
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10 Conclusions 

The performance of muon reconstruction, identification, and triggering in CMS has been stud- 
ied extensively using 40 pb^^ of data collected in pp collisions at = 7TeV at the LHC in 2010. 
These data were used to study several representative muon selections, which were chosen as 
benchmarks covering a wide range of physics analysis needs. 

The distributions of kinematic and identification variables for muons reconstructed from in- 
clusive data samples are generally in good agreement with predictions from simulation over 
the momentum range of pj < 200GeV/c, including variables sensitive to large hit occupancies 
in the muon detectors. Identification efficiencies, for muons with pj larger than a few GeV/ c, 
are above 95% for all selections studied. They are correctly reproduced by the Monte Carlo 
simulations. Misidentification is lower than 1% for the loosest selection and below 0.1% for the 
tightest. In the specific case of several nearby muons, where the tightest selection becomes less 
efficient, an optimized selection was developed to preserve high efficiency with only a minimal 
increase of misidentification. 

Two complementary methods have been developed to evaluate muon momentum scale biases 
and resolution in the range 20 < px < 100 GeV/c, where the momentum measurement is pro- 
vided by the tracker. The average bias A(pt)/Pt in the muon momentum scale was measured 
with a precision of better than 0.2% and was found to be consistent with zero. The relative px 
resolution is between 1.3% to 2.0% for muons in the barrel and better than 6% in the endcaps, in 
good agreement with simulation. The transverse momentum resolution of tracks reconstructed 
using information from only muon detectors is better than 10% in the barrel region. At high 
momenta, the best measurement of muon pj is obtained by selective use of information from 
the muon system in addition to that from the inner tracker. The pj resolution, evaluated in the 
barrel region using muons from cosmic rays, is better than 10% up to 1 TeV/c. 

Algorithms to identify cosmic and beam-halo backgrounds among collision events were devel- 
oped and successfully used in physics analyses of 2010 data. The performance of various muon 
isolation algorithms was shown to be reasonably well modelled by the simulation. 

The muon trigger efficiency for isolated muons is better than 90% over the full t] range, and is 
typically substantially better. 

In this document we have shown that the performance specifications set out for the measure- 
ment of muons in CMS have largely been met. The good performance and detailed under- 
standing of the muon reconstruction, identification, and triggering provides the necessary con- 
fidence in all elements of the chain from muon detection to muon analysis, which is essential 
for searches for physics beyond the Standard Model as well as accurate Standard Model mea- 
surements. 
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